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ABSTRACT 

The r-process nucleosynthesis in core-collapse supernovae (CC-SNe) is studied, with a focus on the 
explosion scenario induced by rotation and strong magnetic fields. Nucleosynthesis calculations are 
conducted based on magneto-hydrodynamical explosion models with a wide range of parameters for 
initial rotation and magnetic fields. The explosion models are classified in two different types: i.e., 
prompt-magnetic-jet and delayed-magnetic-jet, for which the magnetic fields of proto-neutron stars 
(PNSs) during collapse and the core-bounce are strong and comparatively moderate, respectively. 
Following the hydrodynamical trajectories of each explosion model, we confirmed that r-processes 
successfully occur in the prompt-magnetic-jets, which produce heavy nuclei including actinides. On 
the other hand, the r-process in the delayed-magnetic-jet is suppressed, which synthesizes only nuclei 
up to the second peak [A ~ 130). Thus, the r-process in the delayed-magnetic-jets could explain only 
“weak r-process” patterns observed in metal-poor stars rather than the “main r-process”, represented 
by the solar abundances. Our results imply that core-collapse supernovae are possible astronomical 
sources of heavy r-process elements if their magnetic fields are strong enough, while weaker magnetic 
explosions may produce “weak r-process” patterns {A < 130). We show the potential importance 
and necessity of magneto-rotational supernovae for explaining the galactic chemical evolution, as well 
as abundances of r-process enhanced metal-poor stars. We also examine the effects of the remaining 
uncertainties in the nature of PNSs due to weak interactions that determine the final neutron-richness 
of ejecta. Additionally, we briefly discuss radioactive isotope yields in primary jets (e.g., ®®Ni), with 
relation to several optical observation of SNe and relevant high-energy astronomical phenomena. 
Subject headings: nuclear reactions, nucleosynthesis, abundances - supernovae: general - stars: neu¬ 
tron - magnetohydrodynamics (MHD) - neutrinos 


1. INTRODUGTION 

The r-process—the nucleosynthesis process of rapid 
neutron capture compared to /3-decay—is considered to 
be a main astronomical source of nuclei heavier than 
iron, including rare-Earth elements and actinides. As¬ 
trophysical sites for the r-process, which provide suitable 
neutron-rich ejecta in explosive conditions, are still un¬ 
determined, although a large amount of effort by nuclear 
physicists, astrophysicists, and astronomers has been de¬ 
voted since the d awn of modern science (for the most 
recen t status, see, |Arnould et al.|2007[ [Thielemann et al.| 


20111. The lack of knowledge about the astronomical 


sources of r-process elements also causes our incomplete 
understanding of the galactic chemical evolution of r- 
process elements and observed abundance patterns of 
very metal-poor stars. 

Core-collapse supernovae (CC-SNe), which are a main 


Hashimotol 

19951 IWoosIev & Weaver||1995l |ThieIemann| 

et ai. II 1996 

Rauscher et al.|2002 Nakamura et al.||2014p. 


have also been expected to provide reliable conditions 
for the r-process, because the central remnant object, a 
neutron star (NS), consists of extreme neutron-rich ma¬ 
terial. In this context, the prompt explosion mechanism, 
which promptly ejects very neutr on-rich matter, was ex¬ 
amined in early studie s (see, e.g., Hillebrandt et al.|1976 


itself has been completely excluded from possible explo¬ 
sion mechanisms by more sophisticated hydrod ynamics 
simulations w i th detailed microphysics inputs (Rampp] 
' fc Janka]|2000[ Liebenddrfer et al. 12001 ) that employ the 
Boltzmann neutrino transport and realistic nuclear eq 


■equa- 


I transport 
tions of states (EoSs). 

The current standard scenario for CC-SNe (the so- 
called “delayed explosion”) is based on a neutrino¬ 
heating mechanism triggered by convection and hy drody¬ 


namical instabil ities, reviving a stalled shockwave (Bethe 


& Wilson 1985). Following this scenario, the ejecta ap- 
pears to be slightly neutron-rich or even proton-rich in 
a late phase d ue to neutrino absorption (for recent re¬ 
views see, e.g., |Janka||2012| Kotake et al.||201^ |Burro\^ 


2013), which provides insufficient conditions for produc¬ 
ing heavy r-process elements. Recent studies have re¬ 
ported that only some specific cases can lead to successful 
explosions avoiding the strong influence of the neutrino 
absorption. One candidate is the explosion of lighter 
massive progenitors (< 9 M(t,), called the 0-Ne-Mg or 
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mostly produce lighter r-process elements {A < 130). 

On the other hand, the collapsed core of a massive 
star, a nascent proto-neutron star (PNS), has a wind-like 
outflow from the surface driven by strong neutrino emis¬ 
sion. This outflow is called the PNS wind (or neutrino- 
driven wind) and has been regarded as the most reli¬ 
able astrophysical site for the r-proc ess, since the first 
proposal and pioneering studies (e.g., Meyer et al.||1992 


Woosley fc HoffmanI 19921 'I'akahashi et al. 1994| [witti 


et al.||1994| [Takahashi et al.||1994p. However, aside 


from the long-standing problem of high entropies for a 
successful r-process, the latest hydrodynamics simula¬ 
tions with detailed microphysics have revealed that the 


proton-rich (Fischer et al. 

20101 Hiidepohl et al. 

2010 

Arcones fc Thielemann||2(J13|). 'i'heretore, it appears in- 


cessful r-process (revisited by Wanajo 2013[ taking into 


account the recent progress of PNS physics). However, 
more advanced treatments of neutrino reactions in dense 
matter possibly change the electron fraction (Yp) of the 
winds (Martinez-Pinedo et al.|2012 Roberts et al.|2012| 
and optimistically co uld lead to a weak r-process (see, 
also |Wu et ar]|2015[ for the possible influence of neu¬ 
trino osallatiOTsjnn conclusion, the PNS wind scenario 
should be excluded from possible production sites for 
heavy r-process nuclei, including the third park elements 
(A ~ 195) and actinides. 

Although the PNS wind from canonical CC-SNe hardly 
produces heavy r-process nuclei {A > 120), in case of fast 
rotation and strong magnetic fields the physical condi¬ 
tions for the r-process can change. Strongly magnetized 
PNSs, or “proto-magnetars”, have been studied as an al¬ 
ternative scenario of regular (in the spherical symmetry) 
PNS winds (see, e.g., Nagataki fc Kohri|2001 Thompson 


|et al.|20M| [Suzuki &: iNagataki|20U5[|Metzger et al. 

and references therein). Using parameters of rotation 

and magnetic fields based on several PNS wind models, 
more high entropies and/or a rapid expansion of winds 
seems more favorable than conditions in non-rotating, 
unmagnetized P NSs to permit the production of heavy 
r-process nuclei d Vlasov et al.|2014|). Further possible nu¬ 
cleosynthesis signatures of m agnetars (e.g., the produc - 
tion of ®®Ni), are discussed in Suwa & Tominaga|(2015). 

Besides CC-SNe and PNS winds, other important can- 
didates are remnant compact binaries, NS, and NS/black 
hole (BH) pairs, which eject neutron-rich matter by their 
collision. These compact binary mergers (CBMs) have 
also been in vestigated as another possible r-process sce¬ 


nario (e.g. 


Lattimer & Schramm 


et al.||1999a^ . More sophisticated inve stigations (e.^ 


tures and neutrino absorption). Beside this very neutron- 
rich dynamical ejecta of CBMs, wind outflows around 
the central NS and accretion disk are expected to fill 
this underproduction of lighter r-process elements. The 
disk outflow via turbulent and nuclear heating co nsists 
of moderate neutron-rich (1/ > 0.2) matter (e.g. 
zger et ar]|2009a|b [Just et ari|2015( [Perego et al.| [! 


Met- 


Rosswog et al.|2014 1 . Moreover, additional outflows with 
relatively higher driven by neutrino are expected for 
cases of longer durat ion (> 100 ms) of BH format ion from 
a hyper-massive NS ( Metzger fc Fernandez|20M) and the 


accretion disk of BHs (e.g., Wanajo fc Janka||2U12 1. 


The CBM scenario, however, struggles to understand 
the observational signatures in very early galaxies. It 
seems difficult to explain the existence of r-process ele¬ 
ments in the early g alaxy at the extrem ely low metallic- 


ity of [Fe/H] = —3 (Argast et al. 2004), because CBMs 
are delayed due to tne formation of binary NSs and the 
timescale of merging. This implies that, alternatively, 
CC-SNe may contribute galactic r-process elements at 
least in the early universe (at low metallically). There¬ 
fore, potentially multiple astrophysical sites for the r- 
process are required, even if CBMs are the major com¬ 
ponent. Recently, several studies on galactic chemical 
evolution have revisited this problem in different models 


and CC-SNe (e.g., Cescutti fc Chiappini|2014 

Matteuccif 

et al. 

2014 

I'sujimoto fc Shigeyamaj|2014a|b 

|lshimaru| 

et al. 

2015 

van de Voort et al.| 

20151 IShen et al.||2015l 

Wehmeyer et al. 2015). The fu 

1 answer to resolve the 


enhancement history of the early Galaxy requires inho¬ 
mogeneous chemical evolution studies with sufficiently 
high resolution to follow the pollution of the interstel¬ 
lar medium by individual events, which also include all 
possible mixing processes. 

CC-SNe induced by strong magnetic fields and/or fast 
rotation of the stellar core (i.e., magneto-hydrodynamical 
supernovae, MR-SNe)[^ are considered to provide an 
alternative a nd robust astronomical source fo r the r- 
process (e.g., Symbalist^|1984 Cameron 2003). Nucle- 


osynth etic studies were carried out by |]N ishimura et al. 


(|2006|) based on adiabatic MHD simulations, which ex- 
mbit a successful r-process in jet-like explosions. Ad¬ 
ditionally, magnetically driven jets of “collapsar models” 
(BH accretion dis ks) also have been investigated as a site 


of the r-process (Fujimoto et al. 

2006 

Nagataki et al. 

2007[ Harikae et al. 2009 

). Strong^ magneto-rotational 


1974[ Freiburghaus r-pro cess nuclei (Fujimoto et al. 2007 2008 Ono et al. 


Coriely et af.|20lTj|Korobkin et al.|2012 [Bauswein et al.| 

shi 


. . 1 ^ 

2013p , have sh own that the ejecta of CBMs remains ex- 

tremely neutron-rich, producing heavier nuclei than the 
second r-process peak due to strong neutron capture with 
fission recycling. On the other hand, r-process elements 
up to the second peak {A ~ 130) in the dynamical ejecta 
are severely underproduced, because fission fragment dis¬ 
tributions do not extend to isotopes much below this 
mas s number. The most recent studies on NS-NS merg¬ 
ers (Wanajo et al. 2014 Just et al. 2015) address this 
open question, based on the shock-heating ejection sce¬ 
nario with stronger weak interactions (i.e., position cap¬ 


2012 I with a very simple treatment of BH formation. 
Here, we should note that the above noted studies for 
MR-SNe and collapsar models were performed in the ax- 
isymmetric MHD simulations with simplified treatments 
of neutrino transport. 

One important question is whether the earlier re¬ 
sults assuming axis-symmetry also hold in full three- 

® Here, this scenario is distinguished from neutrino- and mag¬ 
netic driven winds from the proto-magnetar, although they may 
have some association (i.e., the jet-like explosion of MR-SNe has 
been considered to be a birth event of magnetars). 

® The disk wind of colla psar m odels were also investigated (e.g., 
[Surman et al.||2006| [2008112014| | as sites of heavy element nucle- 
osynttiesis, mostly based on semi-analytic models. 
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dimensional (3D) simulations (i.e., lead to the ejection 
of jets along the polar axis). MHD simulations in 3D 
with an im proved treatment of ne utrino physics were per¬ 
formed by Winteler et al. (2012) for a ISMq progenitor 
utilizing the initial dipole magnetic field of 5 x 10^^ G 
and the ratio of magnetic to gravitational binding energy, 
£'mag/|bk| = 2.63 X 10“®. These calculations supported 
and confirmed the ejection of polar jets in 3D, attain¬ 
ing ^ 5 X 10^® G and ifmag/|f^| = 3-02 x 10“^ at the 
core-bounce, with a successful r-process that produces 
up to and beyond the third r-process peak {A ~ 195). 
More recent 3D-M HD simulations i n the general rela¬ 
tivistic framework (Mosta et al. ]M4l ), involving a 25Mq 
progenitor with initial magnetic fields of 10^^ G, initially 
lead to jet formation, which afterward experiences a kink 
instability and deforms the jet-like feature. 

Possibly the difference between the two latter investi¬ 
gations in 3D hydrodynamics marks a transition due to 
passing critical limits in the stellar mass and the initial 
rotation and magnetic fields between a clear jet-like ex¬ 
plosion and a deformed explosion. Nevertheless, when 
we see the time evolution of the expl osions with the hy - 
drodynamic instability (see, results in Mosta et al. IMlj ), 
the region with the highest magnetic pressure contains 
essentially the matter corresponding to the initially form¬ 
ing jets before the deformation, which is expected to keep 
the nucleosynthetic features. For this reason, we justify 
continuing with the axis-symmetric simulations instead 
of 3D-MHD in the current study, even if the initial jet 
features will not remain until the completion of the ex¬ 
plosion and the production of r-process elements. We ex¬ 
pect that the overall nucleosynthesis composition is close 
to reality, which is determined by the neutron-richness of 
the ejecta rather than the property of the shock propa¬ 
gation and the timescale of expansion in the outer layers. 

The present study focuses on a remaining problem in 
the MR-SN scenario (i.e., the detailed production and 
ejection mechanism of r-process nuclei). A series of long¬ 
term explosion simu lations, based on the special rela- 


tivistic (SR) MHD (Takiwaki et al.|j2009 Takiwaki & 
Kotake 2011), is adopted tor the nucleosynthesis exam¬ 
ination. i'hese simulations follow the amplification of 
magnetic fields due to differential rotation (winding of 
magnetic fields) and the launch of the jet-like explosion. 
Our MR-SN models include a realistic nuclear EoS, neu¬ 
trino emission from the PNS, and a simplified treatment 
of neutrino transport b y the multi-flavor leaka ge sc heme. 


Following the studies o f Takiwaki et al. (2009) and Taki¬ 
waki & Kotake (2011), we classify our explosion niodels 


into two categories - prompt-magnetic-jet and delayed- 
magnetic-jet explosions - determined by the strength of 
the magnetic field just after core-bounce. We will show 
that this classification is also applicable for investigating 
the production mechanism of r-process nuclei. 

This paper has the following structure. Section 
briefly summarizes basic physics and numerical methods 
for hydrodynamics and nucleosynthesis simulations used 
in the present study. The scenario of MR-SNe and the 
physical properties of explosion models are explained in 
Section The following Section describes the results 
of nucleosynthesis, including the r-process, and we also 
discuss the impacts of the physical uncertainty of MR- 
SNe on nucleosynthesis in detail. Finally, Section is 


devoted to summary and conclusions. 

2. METHODS 

We summarize the basic physics and numeri cal meth¬ 
ods of our M R-SN models, which are based on |Takiwaki| 


et al. (|2009|). Additionally, we describe the TPM, the 
treatment of weak interactions at the nuclear statistical 
equilibrium (NSE) and relevant physical uncertainties, 
and the nuclear reaction network. 

2.1. Pre-collapse Models 

We prepare pre-collapse models with rotation and mag¬ 
netic fields used in the following hydrodynamical explo¬ 
sion simulations. Adopting a 25Mq pre-collapse model 
evolved in the spherical symmetry (ID), we assume sim¬ 
ple analytic formula for the distribution of stellar rota¬ 
tion and magnetic fields at the gravitational collapse, be¬ 
cause the structure and dynamics of the stellar interior 
in the la te phases of evolution are still poorly known 
(see, e.g.. Spruit [2002 for rotation and magnetic fields). 
The key physical properties of our pre-collapse models 
are summarized as follows: 

1. Hydrodynamic structure —A massive rotating stel¬ 
lar evolution model i n the spherical sym metry (E25 
model; described in Heger et al. 2000) is adopted 
as the standard progenitor in this study. This is an 
evolution model that has a 25 Mq zero-age main- 
sequence mass and becomes a Wolf~Rayet star dur¬ 
ing the core He-burning. The thermodynamical 
quantities of the model at the end of evolution (the 
pre-collapse stage; i.e., density, internal energy, en¬ 
tropy, and Ye) are used as initial conditions for hy¬ 
drodynamical simulations. The enclosed mass and 
radius of the “iron core,” which is composed of iron- 
group nuclei (e.g., Fe and Ni isotopes), are 1.69 Mq 
and 2188 km, respectively. 

2. Rotation — As commonly used in CC-SN simula¬ 
tions (see, e.g., Kotake et al.||2003 Takiwaki et al. 


2004), the cylindrical rotation prohle is employed 
tor tne initial rotation, which is expressed by the 
angular velocity in the meridian plane: 


n{x,z) = Do 


An 


^n 


A2 + Ao" ^4 + 


( 1 ) 


where A and Z denote distances from the rota¬ 
tional axis and equatorial plane. The constant 
value of Do represents the strength of rotation, 
while Ao and Zq determine the scale for the re¬ 
gion of the fast rotating core (strength of differen¬ 
tial rotation) in the direction of equatorial plane 
and rotational axis, respectively. The strength of 
the rotation is also described by /3 = ifrot/|bF|, 
which is the ratio of rotational energy (Arot) to the 
absolute value of the gravitational binding energy 
(|VF|). In the present study, the shape of rota¬ 
tion is fixed by the constants of Aq = 100 km and 
Zq = 1000 km, whereas the strength of rotation 
varies as Dq = 151, 76, and 38 in rad s”^^ which 
correspond to /3 = 0.25 x 10“^, 1.00 x 10“^, and 
4.00 X 10“^, respectively. 

3. Magnetic fields — Same as rotation, an analytic 
formula of the initial magnetic fields is employed 
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TABLE 1 

Parameters of Pre-collapse Models 


Name 

Bo 

(G) 

Emag/|1T| 

Oo 

(rad s~^) 

0 

Bll/30.25 

10“ 

2.5 X 10“® 

38 

0.25 X 10-2 

Bll/31.00 

10“ 

2.5 X 10“® 

76 

1.00 X 10-2 

B12/30.25 

1012 

2.5 X 10-1 

38 

0.25 X 10-2 

B12/31.00 

10l2 

2.5 X 10-1 

76 

1.10 X 10-2 

B12/34.00 

10l2 

2.5 X 10-1 

151 

4.00 X 10-2 


Note. — Adopted parameters in equations Q and where 
Bq is the strength of magnetic fields in G, Emag/lW| is the ratio of 
magnetic energy to \ W\, Hq is the angular velocity in rad s“^, and 
/3 is the ratio of rotational energy to \W\. Other scale parameters 
are assumed as constant: Xq = 100 km, Zq = 1000 km, and 
tq = 2000 km, respectively. 

for the simplicity. We assume that the initial mag¬ 
netic fields have a poloidal dominant structure in 
the iron core, which is uniform and parallel to the 
rotation axis inside the core and has a dipole shape 
outside the core. The shape of the magnetic fields 
are expressed by an effective vector potential A in 
the spherical coordinate (r,d,</)): 


Ar(r,9) = Ag{r,9) = 0 , 

A^{r,9) = ^ i-sin^ , 

2 


( 2 ) 


which are defined on the meridian (r, 0)-plane. We 
choose the radial scale tq — 2000 km, which ap¬ 
proximately corresponds to the boundary of the 
iron core. Additionally, we use two different mag¬ 
netic field strengths Bq = lx 10^^ G and 1 x 10^^ G, 
which correspond to i?mag/|kb| = 2.5 x 10“® and 
2.5 X 10“^, respectively. Here, i?mag/|kb| denotes 
the ratio of magnetic energy (ifmag) to the absolute 
value of the total gravitational binding energy. 

All of the adopted parameters in Equations Q and ([^ 
are summarized in Table As shown in this table, we 
arrange five different initial conditions for core-collapse 
and explosion simulations. In the present paper, each 
model is named after its initial quantities of Bq in G and 
j3 = Erot/\W\ (of 10-2), (i.e., “Bll/30.25,” “Bll/31.00,” 
“B12,50.25,” “B12,dl.00,” and “B12/34.00”). The ex¬ 
plosion scenario of MR-SNe and properties of explosion 
models based on the ab ove initial conditions are de¬ 
scribed in Section [3T] and respectively. 


2.2. Special Relativistic Magnetohydrodynamics 

The nucleosynthesis calculations of the r-process, in 
general, require long-term hydrodynamical evolution to 
determine unbound matter (neutron-rich ejecta). For 
this reason, we employ an SR-MHD code, which has 
higher computational efficiency and physical validity of 
dynamics around the speed of light. Previous investiga¬ 
tions based on Newtonian MHD have reported that the 
Alfven velocity unphysically exceeds the speed of light 
when the jet-hea d of the shock outflow has strong mag¬ 
netic fields (e.g., Takiwaki et al. 20041. The Alfven ve¬ 
locity is simply estimated by 


B 


VA = 


y/Mrp 


10^° cm s"^ 


B/1013 G 


a/p/ 10® g cm“3 


(3) 


where p and B are the density and the strength of 
magnetic field, respectively. If we take typical values 
around the jet-head of outflow (i.e., B ^ 10^^ G and 
p ^ 10® g cm”®), the second term reaches unity, where 
va ^ 10^° cm s”^ exceeds the speed of light. In addi¬ 
tion to violating the relativistic causality condition, this 
causes the problem of numerical efficiency, because the 
maximum va limits the stable numerical interval of each 
time step by the Courant-Friedrichs-Lewy condition. 

For this reason, we use an SR-MHD code that has 


been developed for GG-SN simu lations (Takiwaki et al. 


2009 Takiwaki & Kotake 2011), based on an F uferian 


hnite-ditterence method (Stone & Norman 1992). This 
SR-MHD code follows the relativistic torinalisni of IDel 


Villiers et al. (2003), where states of relativistic fluid 
are expressed by hydrodynamic quantities (i.e., density, 
p; internal energy, e; velocity, u®; and pressure, p) at 
each point in the space-time geometry. This code ap¬ 
plies a constrained transport scheme and the method 
of characteristics to manage the induction equation and 
the divergence-free condition, respectively. The Poisson 
equation of self gravity is solved by the modified incom¬ 
plete Cholesky conjugate gradient method. While solv¬ 


ing MHD e quations, we also employ a nuclear EoS (jShen 
et al. |1998 ) in the bases of the relativistic mean field the¬ 
ory. 'The original EoS table is limited for high densities, 
so that we assume that an ideal gas EoS is valid in low 
densities. More details on th e numerical aspect and com- 
puting tests are available in Takiwaki et al. (2009) and 
references therein. 

We perform hydrodynamical simulations on a two- 
dimensional (2D) computational domain with the as¬ 
sumption of axisymmetry. The computational region of 
the code spreads a quarter of the meridian plane with 
equatorial symmetry. We adopt the 2D mesh in the 
spherical coordinates, which have 300 non-uniform grid 
points in the radial-direction (r) and 60 uniform points 
in the polar direction (9), respectively. The radial grids 
cover from 0.5 to 4000 km in the radial-direction, where 
the innermost grid has the finest interval of 1 km. The 
polar grids constantly spread in the range of 0 = 0 - 7r/2 
rad, of which the minimal grid size corresponds to 25 m. 
Additionally, we assume a reflection boundary condition 
for both the polar axis and equatorial plane. The nu¬ 
merical convergence of hydrodynamical simulations was 
performed using these 2D special grids o n MR-SN simu- 
lations, and is described in Section 5 of|Takiwaki et al.| 
([ 200 ^. 

Nucleosynthesis studies on the r-process of MR-SNe 
(Nishimura et al.||2006 Winteler et al.|2012 ) have shown 
that the Ye of ejecta is influenced by the electron and 
positron capture and neutrino absorption in the inner 
region (< 1000 km) during collapse and explosions. This 
implies that weak interactions including neutrino absorp¬ 
tion are significant, as well as for explosion dynamics 
and formation of the PNS. However, the neutrino trans¬ 
port has never been solved in a fully coupled with multi¬ 
dimensional hydrodynamics systems, even by state-of- 
the-art GG-SN simulations due to computational diffi¬ 
culties. Thus, we adopt an approximate method for the 
neutrino cooling of the PNS by means of the neutrino 


multi -flavour leakage scheme (Rosswog & Liebendorfer 


2003). This scheme takes into account the change of 
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TABLE 2 

Initial distribution of tracer particles 


Label 

A^all 

{Nr,Ng) 

Ar 

(km) 

A6 

(rad) 

m 

(Mq) 

^min 

(Mq) 

"^max 

{Mq) 

Low 

2500 

(50, 50) 

79.8 

1.54 X 10“^ 

8.19 X 10-"^ 

1.22 X 10-* 

2.77 X 10-3 

Medium 

10000 

(100,100) 

39.9 

7.72 X 10-3 

2.05 X 10-"^ 

1.07 X 10-5 

7.00 X 10-"^ 

High 

50000 

(500,100) 

7.98 

7.72 X 10-3 

4.10 X 10-5 

4.34 X 10-^ 

1.40 X lO-'^ 


Note. — TVall is the total number of tracer particles. Nr and Ng are the number of particles located in the radial and 6-direction 
(iVall = Nr X Ng). Their intervals are Ar and A6, respectively, fn is the average, rrimin is the minimum, and rrimax is the maximum masses 
of particles in Mq. 


the Yf. by weak-interaction (i.e., electron and positron 
captures on nucleons), and photo-, p air-, and plasma- 
processes (reaction rates are taken by Takah ashi et al. 
1978 Fuller et "^llOSS |Itoh et ar]|1989 ' 1996 b 


The neutrino-spheres ot all flavors in which neutri¬ 
nos are trapped are consistently defined as the neutrino- 
optically thick region by the leakage scheme; whereas, 
emitted neutrinos from the sphere freely escape outside 
the surface (the neutrino-optically thin region). There¬ 
fore, when coupled with the leakage scheme, our MHD 
simulations ignore the effect of neutrino transport and 
reactions in the neutrino-optically thin region (outside 
the neutrino-sphere), where the Ijs changes due to neu¬ 
trino absorption. Including these effects on nucleosyn¬ 
thesis calculations, we take into account the change of Ye 
outside the neutrino-sphere within a post-process TPM, 
which is described in the following section. 

2.3. Tracer Particles with Evolution of Yg 

The MHD code used in the present study employs 
grid-based (Eulerian) hydrodynamics, of which fluid dy¬ 
namics is given by the time evolution of fluid on spa¬ 
tial numerical mesh. The particle-based (Lagrangian) 
hydrodynamics, in contrast, make use of the initial posi¬ 
tion (or the particle’s number) to describe fluid motion, 
which is more reliable to follow the abundance evolu¬ 
tion of ejecta. We adopt the so-called tracer/test par¬ 
ticle method (TPM), which has been used to convert 


of Lagrangian particle motion (see, e.g., Nagataki et al. 

1997 Travaglio et al. 

2004| Seitenzahl et al. 

20i0|. h'or 


sion of the previous one (used in Nishimura et al. 2006 
Wanajo et al.||2014 ). 


2.3.1. Numerical schemes 

We briefly summarize the numerical methods (i.e., the 
time integration and spacial interpolation) of the TPM 
code. The time integration is based on the two variables 
(2D) predictor-collector method with the second order 
accuracy. At each computational time (t": the nth time 
step), we adopt the bilinear interpolation for converting 
the grid-based physical quantity into tracer particles. In 
polar coordinates (r, 9), the TPM code calculates the po¬ 
sition of the tracer particle at (i.e., (r"+^,0"+^)), 
using the previous position (r", 0"’) and the velocities in 
the r-direction (vr) and 0-direction (vg). This relation is 
expressed by 


r"+i = r" -b u!" At” 


gn+l 


(4) 


where At”(= - t”), , and are the time 

interval, modified-radial-velocity, and modified-angular- 
velocity, respectively. Here, u*" and are functions 
of velocities at two previous time steps and t”, es¬ 
timated by the predictor-collector scheme. Additionally, 
the refraction boundary is assumed for the rotational axis 
and equatorial plane, which is the same as the original 
MHD hydrodynamics simulations. 

For the initial position of particles, we adopt a uni¬ 
form distribution in space at the beginning of the simu¬ 
lation (on the pre-collapse model), with different particle 
masses. Tracer particles are initially located in a quarter 
of the meridian plane, whose domain is r = 0.5-4000 km 
and 0 = 0 - 7r/2 rad. The innermost initial radius of trac¬ 
ers is r = 50 km. We prepare three different sets of tracer 
particles in different spatial resolution for the purpose of 
a numerical convergence test. Tablej^shows their param¬ 
eters and relevant physical quantities. The total number 
of particles for each set is expressed by Aaii = Nr x Ng, 
where W and Ng are the number of particles in the r- 
direction and the 0-direction, respectively. Additionally, 
TO, TOminj and TOmax are the average, maximum, and min¬ 
imum value of masses. The mass of each tracer particle 
is conserved during the motion, of which the initial inner 
particles have larger values compared to ones in outer 
layers. The particle sets are labeled low, medium, and 
high expressing their numerical resolution. The numer¬ 
ical convergence of our TPM estimation is discussed in 
Sections and focusing on the Yg of tracer particles 
and r-process nucleosynthesis, respectively. 


2.3.2. The Evolution of Yg Based on the TPM 

The physical value of a tracer particle basically fol¬ 
lows the evolution of original hydrodynamical model that 
neglects the effect of neutrino absor ptio n outside the 
neutrino-sphere (described in Section 2.2). These neu¬ 
trino reactions, however, may be crucial for determina¬ 
tion of the final Yg, which is important for the resulting 
r-process abundances. Thus, we will include neutrino- 
absorption reactions on the TPM during expansion. For 
this purpose, we define times from the onset of core¬ 
collapse (i.e., tjy and t^se) for each particle, focusing on 
the moment of escape from the neutrino-sphere and NSE 
state. Here, t^, is the moment just before the particle es¬ 
capes from the neutrino-sphere and tnse is the time at 
the end of NSE that the temperature drops to 9 x 10® K 
(= 9 GK). 

Figure provides a basic picture of the ejection pro¬ 
cess for a typical tracer particle with corresponding time 
variables {tn and fnse)- We can safely assume tnse > tv 
for all particles because of the property of our explosion 
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Fig. 1.— A schematic picture of tracer particle motion. On the 
meridian plane, the tracer particle that is initially located outside 
of the core (lower right) is finally ejected along the rotation axis 
(upper left) via the convective motion inside the neutrino-sphere 
(lower left). During the motion, the particle reaches the edge of 
the neutrino-sphere (t = f^) and escapes form the high temperature 
NSE region (i = tnse)- 


models, that is, the temperature in the neutrino-sphere 
is high enough for the NSE, which agrees with general 
features of CC-SNe. Using these two specific times as 
boundaries, the ejection history of a tracer is divided 
into the following three phases: 

1. Phase 1 {t <t^): This is the stage that the central 
core contracts and the surrounding matter falls, 
triggered by the gravitational collapse. The tem¬ 
perature exceeds 1 MeV (~ 11.6 GK), where heavy 
isotopes are destroyed into free nucleons by photo¬ 
dissociations (the NSE state achieves). On the 
other hand, densities in the inner core also increase 
and become high enough to trap neutrinos inside. 
In this phase, we calculate the evolution of tracer 
particles, adopting the original hydrodynamical ex¬ 
plosion models in the weak equilibrium based on 
the Leakage scheme. 

2. Phase 2 < t < tnse)- The tracer particle has 

already escaped from the neutrino-sphere, but the 
temperature is still high enough (T > 9 GK) for 
the NSE. Although hydrodynamic simulations ig¬ 
nore neutrino absorption, we include the effect on 
changing by means of a post-process treatment. 
Using the geometric property of the tracer particle 
and neutrino-sphere, we take into account electron 
and anti-electron neutrino captures on nucleons. 
The reaction rates of neutrino absorption are de¬ 
termined by the luminosity, the mean energy, and 
the radius of neutrino-spheres. The isotopic abun¬ 


dances are obtained by solving the NSE (nuclear 
Saha) equation that is a function of the tempera¬ 
ture and electron density (determined by the mat¬ 
ter density and Yg). 

3. Phase 3 (tnse < t)-. After the temperature of the 
tracer particles drops to 9 GK, the NSE becomes 
invalid, where all relevant nuclear reactions are 
activated. At this stage, we calculate the abun¬ 
dance evolution using a full nuclear reaction net¬ 
work. The network calculation follows all relevant 
nuclear reactions involved in the the r-process and 
other explosive nucleosynthesis, along with thermo¬ 
dynamical evolution of the tracer particle. 

As described, we simply adopt the Y^ from hydro¬ 
dynamics simulations in the neutrino-optically thick re¬ 
gion in the weak equilibrium from the onset of the core¬ 
collapse to t = (Phase 1). For the optically thin region 
(Phases 2 and 3), where tracer particles escaped from the 
neutrino-sphere, we consider the following weak interac¬ 
tions for nucleons: 


p -I- e ^ n -I- r'e , 
n -I- e''" ^ p -I- Ue , 


( 5 ) 


where p, n, i/g, and 14 represent protons, neutrons, elec¬ 
tron neutrinos, and anti-electron neutrinos, respectively. 
While the matter is still in the NSE (the temperature 
exceeds 9 GK), we calculate the time evolution of Yg, 
described by the equation: 

dY 

= ~(^pe- + Ap!7e)T^ + (Ane+ + An;y3)(l — Yg) , (6) 


where Ape-, A„e+, Amy^, and App^ denote forward and 
reverse reaction rates in Equation (™, determined by a 
couple of density and temperature. \Ve adop t reaction 
rates for electron and posi tron captures from |Langanke| 


& Martmez-Pinedo (2001) given in a grid of tempera- 
tures and electron densities. For their inverse reactions, 


i.e., neutrino abso rptions, simple analytic formulae (Qian 


& Woosley 1996) are used as functions of the neutrino- 
sphere radius and the luminosity and mean energy of 
neutrino emission from the PNS. 

The duration of hydrodynamical simulations (^ 
100 ms after the core-bounce) is insufficient to follow 
all relevant nucleosynthetic processes, even for the r- 
process. For the r-process, neutron captures and /3- 
decays continue several hundred seconds (^ 1 - 100 s) 
during expansion. We obviously do not have the time 
evolution of hydrodynamical values based on MHD sim¬ 
ulations in the phase of r-process nucleosynthesis. How¬ 
ever, in such a late stage of the explosion, the dependence 
of the hydrodynamic evolution (the expansion timescale) 
on nucleosynthesis becomes comparably weak, so that 
we can safely assume a very simple expansion law. Ad¬ 
ditionally, the most important variable Yg, which deter¬ 
mines neutron-richness, immediately reaches the asymp¬ 
totic value after the NSEH (Phase 3), although the r- 
process itself changes Yg due to /3-decay in later stages 
{t > 100 ms after the bounce). 


^ We will dis cuss the evolution of Yg and other physical values 
in Section|3.2.2| Figure [ t] shows a typical time evolution of Yg for 
selected trajectories. 
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Fig. 2.— Abundances of heavy elements at NSE. The total mass 
fraction of heavy nuclei (AIa), defined by their mass number, A > 
6, are plotted in a typical density range. The cases of Yg = 0.05, 
0.10, 0.15, and 0.20 are calculated with a constant temperature 
T = 9 GK. 


Here, we use an analytic extrapolation to fill the gap 
between the time duration of the simulation and the need 
for nucleosynthesis simulations. Assuming adiabatic ex¬ 
pansion after the end of the simulations (tf), we have the 
constant velocity v(t) = Vf, radius r(t) = rf + (t — tf) Vf, 
density p(t), and temperature T(t), described by 


P(i) = Pi 
T{t) = T{ 


Tf 

r{t). 

Tj 

r{t) ’ 


n 3 


(7) 


where tf, Vf = v(t[), rf, pf, and Tf are the final velocity, 
radius, density, and temperature of the hydrodynami- 
cal simulation, respectively. Although this extrapola¬ 
tion may not be realistic for modeling the hydrodynamics 
of explosion models, the impacts on nucleosynthesis by 
physical uncertainties are expected to be negligible be¬ 
cause our hydrodynamical simulations follow time evo- 
lutio n at least until the end of NSE, satisfying tf > tnse 
(see, Nishimura et al. 2006|. On the other hand, hy- 
drodynamical instabilities may change the dynamics in 
the early stage of Phase 3 100 ms), which launched 


jets that ar e deformed due to the kink instability (Mosta 


et al. 2014) and experienced thermal energy deposition 


by magnetic reconnection. We quantitatively discuss the 
impacts of the kink instability on the r-process in Sec¬ 
tion [4701 


TABLE 3 

Properties of the neutrino-sphere and emission 


Model 

Rue 

Ri^e 




(km) 

(km) 

(10®2 erg s-i) 

(10®^ erg s“^) 

Bll/30.25 

94.0 

79.5 

4.33 

4.76 

Bll/31.00 

102 

79.6 

4.76 

1.77 

B12/30.25 

98.6 

76.6 

3.87 

3.48 

B12/31.00 

101 

78.4 

3.69 

1.53 

B12/34.00 

104 

80.0 

3.52 

6.78 

Note. 

— The 

radiuses 

of neutrino-sphere; 

and R, 7 ^ for 


the electron neutrino and anti-electron neutrino, respectively. The 
luminosities of neutrino emission; Li,^ and after the neutrino 
burst for the electron neutrino and anti-electron neutrino. 


isotopic abundances by solving the NS E (nuclear Saha) 


equati on with weak interactions (see, Nishimura et al. 


2012b I. 

Hy solving the NSE equation, we obtain isotopic abun¬ 
dances at a given set of the temperature, density, and Yg. 
Figure shows the total mass fraction of heavy elements 
(Xa) in NSE, which is the summation of heavy (A > 6) 
isotopes. We calculated abundances in different values of 
Ye, with a constant temperature, 9 GK, the lower bound 
of the NSE. The chosen range of densities and Yg covers 
typical values of neutron-rich ejecta from MR-SNe, where 
the density p < 10^° g cm“^ and Yg = 0.05,0.10,0.15, 
and 0.20. As shown in the figure, heavy elements are neg¬ 
ligible when densities are below ~ 10® g cm“® for any Yg. 
At high densities with low Yg, the values of Xa are still 
less than ^ 10 %, where p > 10® g cm“® with Yg ^ 0.05 - 
0.10, respectively. Therefore, we can assume that heavy 
elements have minor effects on weak interactions. 

In addition to electron and positron captures on nucle¬ 
ons, neutrino absorption is taken into account outside the 
neutrino-sphere (P hase 2 and 3; in Figure!^ . We use the 
absorption rates of|Qian & Woosley| ( 199^, determined 
by the distance from the surface of a neutrino-sphere 
and the physical properties of neutrino emission from 
the PNS, (i.e., the mean energy and luminosity). These 
relevant quantities are consistent with the hydrodynami¬ 
cal explosion models, which are simulated with a simpli¬ 
fied treatment of the neutrino transport by the Leakage 
scheme. The adopted values of the radius of neutrino- 
sphere and the luminosity of neutrino emission are sum¬ 
marized in Table for all explosion models, whose un¬ 
certainties are discussed in more detail in Section 13.2.31 
Using these reaction rates, Equation ®> is solved to cal¬ 
culate the evolution of Yg for each tracer particle along 
the time evolution of the temperature and density. 


2 . 4 . Nuclear Statistical Equilibrium 

During early stages of the explosion (Phase 1 and 2; 
in Figure [l|), most of ejecta are in the state of NSE 
(r > 9 GK], of which isotopic abundances are immedi¬ 
ately determined by the temperature and electron num¬ 
ber density (the product of the density and Yg). If the 
ejecta escaped from the weak equilibrium region (inside 
the neutrino-sphere), the value of Yg changes due to weak 
interactions despite being in the NSE, in which all nu¬ 
clear reactions balance with their inverse reactions. In 
the present study, we assume that the NSE is valid if 
the temperature exceeds 9 GK for any electron num¬ 
ber densities. Above this temperature, we calculate NSE 


2 . 5 . Nuclear Reaction Networks 

We utilize a full nuclear reaction network, which con¬ 
sists of all relevant isotopes and reactions to the r- 
process, for nucleosynthesis calculations. We use an ex¬ 
tension of the previous code, which h as been described in 
detail (see Nishimura et al.||2012a|b ). This network con¬ 
sists of 4423 isotopes from the proton and neutron up to 
fermium (which has an atomic number Z — 100), which 
includes proton-rich isotopes as well as neutron-rich ones 
far from /3-stability isotopes. This reaction network in¬ 
cludes decay channels, two- and three-body reactions, 
electron and positron captures on nuclei, and screening 
effects, as shown in the following equation for the specific 
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ith isotope, denoted by that is 
dt 


^ =j:k,x,y, 


+ ^ M, fc pNav O', k) Yj Yfc 

j, k 

+ J'. k, I 0, k, 1) Y, Yk Yi 

j, k, I 


( 8 ) 


where Nav is Avogadro’s number, and Afij, 
and are proper accounting numbers avoiding 

overlapped-counting, defined by Afij = Ni, Afij^k = 
A'i/n”’Ti|A‘j„|!, and Afi^j^kj = Nij^\l, respec¬ 
tively. Ni is the factor that denotes the variation of the 
ith isotope in all decays and reactions. Additionally, Xi 
is the decay constant, an d (?, k) and ( 7 , fc, 1) are the cross- 
sections (for details, see, Hix fc Thielemann||1999 ). 

At each time step, our network code solves a large lin¬ 
ear equation system, which is based on Equation 
by matri x inversion u sing a sparse matrix solver, UMF- 
PACKpl (Davis 20041. We use a semi-implicit method 
coupleu with the iNewton-Raphson method for the time 
integration. Each thermonuclear reaction rate is given 
as a function of the temperature and density, where 
the thermodynamical evolution of a tracer particle is 
adopted. We choose experimental-based nuclear masses 
provide d by a series of experim ental nuclear mass eval¬ 
uation (Audi & Wapstra 1995). Additionally, theoret¬ 
ical predictions for nuclear masses, reaction rates, and 
beta-decays are applied, based on the finite-range d roplet 
model (ERDM) mass formula ( Moller et al.||1995 ). 

The majority of reaction and decay rates have been 
taken from the compilation of REACLIB databas^ 
( Rauscher fc Thielemann||2000 ) and its extension, JINW 
REACLH^^*^! (ICyburt et al. 2010), in which theoretical 
ones are mainly based on h'RUiVl. Physical uncertainties 
of , 0 -decay half-lives and the effects of recent experimen¬ 
tal result s on the r-process in our network code were dis¬ 
cussed in Nishimura et al. (2012a). We partially utilized 
experiment al-based reaction rates from NACRE (Angulo 
et al. 1999) for lighter nuclei. Eor heavy elements that 
we take into account spontaneous and /3-delayed fis sion, 
we assume an empirical f ormula of fission fragments (Ko- 
dama fc Takahashi|1975 ). Eurther details of adopted iso- 
topes, nuclear reactions, and fissions have been explained 
in our previ ous papers (Nishimura et al.|2006[ Fujimoto 


et al. 2008). We take into account weak interactions, 
whicli are significant to determine the proton/nucleon ra¬ 
tio Ye, as well as nuclear reaction rates described above. 
For the electron and positron capture by heavy nuclei 


including iron-grou 

p isotopes, we use tabled reactions 

(Fuller et al. 

1980 

1982 

Langanke & Marti'nez-Pinedo 


also taken into account, which is explained in the previ¬ 
ous section. Unlike other nuclear reactions that are de¬ 
termined by a combination of the density, temperature, 
and Ye, neutrino-absorption reactions are also depend on 
the distance from the neutrino-sphere and the properties 
of neutrino emission from the PNS. We adopt these quan- 

® http://faculty.cse.tamu.edu/davis/suitesparse.html 
® http://nucastro.org 

https://groups.nscl.msu.edu/jina/reaclib/db 


tities based on our explosion models, which is the same 
manner of estimating Ye during the NS E sta te outside 
the neutrino-sphere (described in Section 2.4). 


3. MAGNETOROTATIONAL SUPERNOVAE 

We describe the explosion scenario of MR-SNe and the 
ejection mechanism of neutron-rich matter, based on our 
explosion models. We focus, in particular, on the evolu¬ 
tion of Ye that directly relates to the following r-process 
nucleosynthesis. Additionally, we discuss the uncertainty 
of Ye in our calculations due to the physical properties 
of a PNS. 


3.1. Explosion scenario 

The mechanism of magnetorotational induced SNe 
(MR-SNe), associated with a jet-like explosion is clearly 
distinguished from the neutrino-heating mechanism of 
canonical CC-SNe. In principle, the dynamics of fast ro¬ 
tating iron cores has common properties from the onset of 
gravitational collapse until the core-bounceF^even if the 


core is significantly magnetized (see, Kotake et al.| 2003 
and references therein). The strength ot the initial rota¬ 
tion in an iron core, which is determined by /3 = T/|IU|, 
governs dynamics before the core-bou nce, where the an¬ 
gular momentum transfer is ignorable (Yamada & Sawaij 
2004). Strong magnetic fields in the PJvS become signif¬ 
icant for explosion dynamics after the core-bounce. 

The dynamics of MR-SNe, which is driven by rapid ro¬ 
tation and strong magnetic fields, may have two different 
types, characterized by the amplification process of the 
magnetic fields around a PNS. Following the study of 
Takiwaki et al. (2009), we classify our explosion mod- 
els into two categories: the prompt-magnetic-jet and 
delayed-magnetic-jet explosions, or simply the prompt-jet 
and delayed-jet, respectively. The former type, prompt- 
jets, promptly explodes just after the core-bounce asso¬ 
ciated with an energetic bipolar-jet explosion, due to a 
strong magnetic pressure by rapid rotation and strong 
magnetic fields at the pre-bounce stage; the latter case, 
delayed-jets, is a relatively weak explosion with a more 
collimated jet-like explosion. The delayed-jets have mod¬ 
erate initial rotation and magnetic fields compared to the 
prompt-jets, and long-term duration until the launch of 
a jet-like outflow after the core-bounce. In the following 
we summarize the characteristic properties of these two 
types for MR-SN models; pictures of 3D hydrodynamical 
structure (in axisymmetry) are shown in Figure]^ 

Prompt-magnetic-jets. —The pre-collapse core has 
strong initial magnetic fields with rapid rotation, and the 
wound-up process during the collapse enhances the mag¬ 
netic fields around the central regionj^ Strong toroidal 
magnetic fields above 10^® G around tEe surface of PNS 
launch a strong shockwave that overcomes the ram pres¬ 
sure of falling matter. This jet-like shockwave runs along 
the rotational axis, so that the outer layers finally ex- 
pload in the polar direction; whereas, shock propaga¬ 
tion in the equatorial plane immediately diminishes and 


The gravitational collapse of the iron core (formation of the 
PNS) is triggered by the same physical processes as standard CC- 
SNe, which are the photo-dissociation and electron capture. 

The “plasma-/?” (this is not /? = T/\W\), which is the ratio 
of the matter pressure to the magne tic pressure, decreases to 10“^ 
in this region ([Takiwaki et al.|2009[l. 
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Fig. 3.— The structure of MR-SNe with a bipolar-jet explosion around the PNS. The entropy per baryon (3-15 fee baryon“^) and 
magnetic field lines (white lines) are drawn in the spatial range of 700 km. The left and right panels correspond to B12/31.00 (prompt-jet 
model) and Bll/30.25 (delayed-jet model), respectively. 


TABLE 4 

Physical properties of explosion models 


Name 

Category 

^bnc 

(ms) 

^del 

(ms) 

^exp 

(ms) 

^fin 

(ms) 

Iljet 

(C) 

-£^exp 

(10^® ergs) 

-^ej,r 

(10-3 Mq) 

(^,r) 

Bll/30.25 

delayed 

244 

48 

72 

103 

0.12 

0.05 

2.68 

0.281 

Bll/31.00 

prompt 

235 

10 

27 

58 

0.27 

0.23 

1.62 

0.192 

B12/30.25 

prompt 

244 

0 

32 

57 

0.20 

1.3 

3.69 

0.186 

B12/31.00 

prompt 

235 

0 

20 

77 

0.27 

1.4 

4.05 

0.221 

B12/34.00 

prompt 

292 

0 

25 

33 

0.20 

1.0 

5.32 

0.181 


Note. — func is the duration of the collapse, which is from the onset of core-collapse until the core-bounce. The other three timescale 
are measured from the core-bounce (tbnc)- fdel i® the raunch of outward shock, tgxp is the shock reaching the distance of 1000 km from the 
PNS, and is the end of simulation. The velocity of the shock at texp (1000 km from the PNS) is Ujeti with a corresponding explosion 
energy of Eexp- ^ej,r 0^e,r) Eire the mass and the average Yg of very neutron-rich ejecta (Y). < 0.3), respectively. 


fails to explode. Only an energetic bipolar-jet explosion, 
which is collimated by strong magnetic fields (the left 
panel of Figure [^, propagates. As shown in the figure, 
the maximum i^ue of entropies inside the jet reaches 
^ 15 fee baryon”^; the central region has lower values 
5 fee baryon”^). 

Delayed-magnetic-jets .—This explosion mechanism is 
caused by comparatively weak initial rotation and mag¬ 
netic helds, which takes longer duration from the core¬ 
bounce for a successful explosion. After the core-bounce, 
the outward shockwave in the direction of the rotational 
axis stalls, as does the the shockwave in the equatorial 
plane, because of a weaker magnetic pressure around a 
PNS. This stalled shock in all directions begins to os¬ 
cillate until ^ 10-100 ms after the bounce. Although 
this oscillation is diminishing and reaches a nearly hy¬ 
drostatic state, the magnetic fields behind the shock 
grow due to magnetic helds wrapping, where the plasma- 
P decreases (i.e, the magnetic pressure becomes pre¬ 
dominant). When the toroidal component of magnetic 
helds exceeds 10^® G, the stalled shockwave revives and 
propagates along the rotational axis. Comparing to the 


prompt-jets, the shape of the delayed-jets is more colli¬ 
mated but less energetic. As shown in the right panel 
of Figure magnetic held lines are strongly wrapped 
around thePNS. 

In the present study, we investigate nucleosynthetic 
properties of the MR-SN, mainly focusing on the mech¬ 
anism of explosion and the ejection process of neutron- 
rich material. Investigations of the MR-SN are also an 
important subject of numerous unsolved problems (i.e., 
macroscopic and microscopic amplihcation processes of 
magnetic helds, neutrino bursts, gravitational wave de¬ 
tection, energetic supernova events, and the study of 
GRBs). The predicted properties of gravitational wave 
signals and neutrino bursts are appar ently diherent be¬ 


tween the prompt- and delayed-jets (Takiwaki & Ko- 


take 2011). Additionally, the effects ot anisotropic neu¬ 


trino bursts and inhuences to the neutrino-oscil lation via 
the jet-like e xplosion have been investigated by|Kawagoe| 


et al. (2009). For more details of those several studies 
and relevant issues to MR-SNe, see Kotake et al. (2012) 
and references therein. 
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Fig. 4. — The dynamics of explosion models expressed by tracer particles for B12/94.00 (left) and Bll/30.25 (right). The distribution of 
particles is shown in a quarter of each panel at the moment the jet-head reaches 1000, 2000, (SOOO, and 4000 km, respectively. Corresponding 
times after the bounce are shown in each panel: which are 10, 25, 52, and 100 ms for prompt-jets and 50, 66, 84, and 109 ms for delayed-jets, 
respectively. The color scale of each particles indicates Ye,nse, Ye at the end of NSE. 


3.2. Explosion Models 

Explosion models of the MR-SN are simulated em ploy- 
ing the grid-based Eulerian code described in Section [2^ 
On the other hand, the particle-based Lagrangian hydro- 
dynamic evolution b ased on the post-process scheme is 
utilized (see, Section 2.31 to follow the abundance evolu¬ 
tion of the models. Several key physical quantities of the 
explosion models are summarized in Table of which 
names cor resp ond to the pre-collapse models described 
in Section |2.1[ The physical properties of dynamics are 


ejecta, {Ye), is defined by 


basically based on the original hydrodynamics models, 
although we adopted the TPM for estimating ejecta. 

The duration of collapse, tbnc, is a time interval in ms 
from the beginning of collapse to the core-bounce. The 
other three timescales are measured from thnc^ ^dei is the 
launch of the MHD jet (delay time), <exp is the jet reach¬ 
ing the distance of 1000 km from the PNS, and tfin is 
the end of simulation. Focusing on the outward jet, the 
velocity is Uexp and the corresponding explosion energy 
is Eexp at the time t = tgxp- The total amount of ejected 
matter, estimated by the post-process TPM, is denoted 
by Mej. There is an alternative mass, denoted by Mej^r, 
which is the total mass of neutron-rich ejecta (producing 
heavy r-process elements) for Ye < 0.3. Finally, the av¬ 
erage electron fraction only for these very neutron-rich 


all 

(Ye, r) = ^ rUi Ye, i (for Ye < 0.3) (9) 


where rrii and Ye, i are the mass and electron fraction 
(at the end of NSE) of individual tracer particles, re¬ 
spectively. This is an index expressing the strength of 
the r-process that lower values produce heavier nuclei, 
mainly because the matter of y(, < 0.3 contributes to the 
production of heavy r-process el ements in low entropy 


Nishimura et al. 2006 


conditions of MR-SNe (see, e.g, _ 

Fujimoto et al. 1120081 INishimura et al.||2012b 
et'kl.|2012] . - -- 


W inteler 


As described in the previous section, prompt-magnetic- 
jets, in principle, have a shorter dynamical timescale, ex¬ 
pressed by tdei and texp, wit h larger ript, and Efixn than 


delayed-jets (for details, see, Takiwaki et al. [[2009 1. Ad¬ 
ditionally, some physical quantities in Table provide a 
clue to understanding not only the dynamics of the ex¬ 
plosion, but also resulting r-process nucleosynthesis. Fo¬ 
cusing on tdei, the delayed-jet model has a larger value 
than prompt-jet models until the launch of the shock. 
Flere, we can predict a stronger progress for the r-process 
for prompt-magnetic models than the delayed-jet model 
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Fig. 5.— The track for ejected tracer particles and the path of tracer particles for a delayed-jet model. Bll/30.25 (upper) and two 
prompt-models, B12/31.00 (medium) and B12/34.00 (lower), are illustrated in different spacial scales. Each particle is labeled by Ye at the 
end of NSE. The central filled (yellow) circular area describes the neutrino-sphere. 


by comparing with although we need precise net¬ 
work calculations to determine the detailed nucleosyn¬ 
thesis yields. 

3 . 2 . 1 . Dynamics of Tracer Particles 

The time evolution of tracer particles for the ejected 
matter is shown for the both of the prompt-jet and 
delayed-jet in Figure The Fg^nse, which indicates Fg 
at the end of NSE, is mso shown for each particle. Each 


panel shows the evolution of the explosion model in dif¬ 
ferent times, for B12/14.00 (left) and Bll/30.25 (right), 
representing the prompt-jet and delayed-jet, respectively. 
A quarter of each panel corresponds to the distribution 
of tracer particles at a given time, in the time sequence, 
which goes through upper left, upper right, lower right, 
and lower left, respectively, evolved in clockwise. The 
time after the bounce is shown in each panel: 10, 25, 
52, and 100 ms for the prompt-jet and 50, 66, 84, and 
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Fig. 6.— The initial position of tracer particles on the meridian 
plane. Ejected particles are plotted with Te, nse at the end of 
NSE. We adopt one delayed-jet (Bll/30.25) and two prompt-jets 
(Bll/31.00 and B12/31.00). 


109 ms for the delayed-jet, respectively. The value of 
Ye distributes in the range of 0.1-0.4 and 0.2-0.4 for the 
prompt model and delayed model, respectively. 

As clearly seen in Figure |4j the prompt-jet tends to 
have a broader structure, of which the jet-head reaches 
^ 3000 km at 25 ms after the bounce, whereas the colli¬ 
mated jet of the delayed model takes ~ 100 ms to reach 
3000 km. Focusing on the Yg, we can see that neutron- 
rich matter (i.e., Yg < 0.3), is ejected in jets along the ro¬ 


tational axis for all the models. Neutron-rich tracer par¬ 
ticles in the prompt-jet model are expected in the body 
of the jet-like outflow, and are continuously launched fol¬ 
lowing the initially ejected matter. The delayed model 
has a different feature: which moderate neutron-rich par¬ 
ticles are mostly populated in a jet-head region, and are 
ejected by the primary jet-like outflow. By only focusing 
on bulk motion of tracer particles, however, we cannot 
fully understand the ejection mechanism of neutron-rich 
matter. 

Therefore, we describe detailed ejection processes us¬ 
ing selected tracer particles with a path on the merid¬ 
ian (X-Z) plane, as illustrated in Figure Three dif¬ 
ferent explosion models are adopted: the delayed-jet 
model, B11/30.25, and two prompt-jets, B12/31.00 and 
B12/34.00. Each explosion model has three panels in 
different spacial scales, and a central area with a ra¬ 
dius of ^ 100 km corresponding to the surface of the 
neutrino-sphere. Each figure has paths of six selected 
tracer particles that follow the history of movement la¬ 
beled by Ee.nsej be., which is the constant of Yg at the end 
of NSE (T = 9 GK). All presented trajectories, in prin¬ 
ciple, fall toward the center during core-collapse and es¬ 
cape along the rotational axis while experiencing chaotic 
convectional motion around the central core. 

The type of prompt-jets represented by B12/31.00 and 
B12/34.00 has an energetic jet-like outflow with a wider 
opening angle, as seen in Figure]^ Tracer particles from 
these explosion models may have low Yg^-g^ along the 
rotational (Z) axis. As shown in Figure M the most 
neutron-rich ejecta in a range of Fe.nse = - 0.2 es¬ 

cape after they deeply fall into the PNS. For the case 
of B12/34.00, which is the most energetic model, a large 
amount of particles are ejected from the central region 
10 km) inside the PNS. By the comparison with 
B12/31.00 and B12/34.00, neutron-rich ejecta tend to 
move the central region, which requires a higher explo¬ 
sion energy to eject them. On the other hand, moderate 
neutron-rich (T^,nse = 0.2-0.3) ejecta are predominant in 
the delayed-magnetic-jet explosion (B11/30.25). For the 
delayed model, trajectories fall into around the PNS core 
and then escape from the core along the rotational axis 
after a long duration 50 ms) of convective motion. 
This duration corresponds to the enhancement process 
of magnetic fiel ds by field wrapping d ue to strong differ¬ 
ential rotation (Takiwaki et al. |2009 ), which results in a 
higher magnetic pressure to overcome the ram pressure 
of falling matter. 

All ejected particles of our explosion models finally es¬ 
cape along the rotational axis, although the convective 
motion in the central area is quite complicated. In order 
to see the initial position of neutron-rich ejecta, we illus¬ 
trate the distribution of tracer particles at the beginni^ 
of collapse with Fg.nse in Figure 1^ Similar to Figure 
we selected Bll/30.25, Bll/31 . OChand B12/31.00, repre¬ 
senting the delayed-magnetic jets and prompt-magnetic 
jets, respectively. The radius of the iron core for the 
progenitor model is ~ 2000 km. Thus, this figure shows 
that ejected particles, in which is significantly be¬ 
low 0.5, are initiated in the iron core. Particles ini¬ 
tially located outside the core along the rotational axis 
are ejected, avoiding the progress of weak interactions 
(i.e., electron-captures), which remain their initial Yg 
(i.e., b^^nse ~ 0.5). For the delayed model, the low- 
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Fig. 7.— Time evolution of the Ye for the selected tracer particles from B12/31.00 (prompt-jet) and Bll/30.25 (delayed-jet). Each line 
begins at the onset of collapse, of which time is measured from at the core-bounce. Thick lines are results taking into account all weak 
reactions between the nucleons and neutrinos, whereas thin lines are results ignoring neutrino absorptions. 


B12P1.00 (prompt) 



- 0.1 0 0.1 
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time after bounce, s 
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Fig. 8.— Time evolution of hydrodynamical values for B12/31.00 (prompt; left) and Bll/30.25 (delayed; right). Each panel includes the 
temperature in K (top), density in g cm“® (medium), and radius in km (low) for selected tracer particles. The condition for the NSE 
(T ~ 9 GK) and the boundary of the neutrino-sphere are highlighted. 
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Fig. 9.— The maximum density and temperature of ejecta with 
Ye for a prompt model, Bll/31.00 (upper), and a delayed model, 
Bll^O. 25 (lower). The ye,nse of each particle is shown in the given 
color range. 


est l^^nse 0.3, which was initially located 1000 km 
from the center, spreads in the direction of 45° from 
the axis. By contrast, prompt-jets have an initial dis¬ 
tribution that from the axis to the equatorial direction, 
especially neutron-rich ejecta that is Ve.nse 0.1, which 
initially locates not only the Z-direction, but also the X- 
direction near the equatorial plane. For all of the models, 
tracer particles initiated within the 500 km of the iron 
core immediately fall into the center and finally construct 
a PNS, which has never been ejected. 

3 . 2 . 2 . Ejection Process of Neutron-rich Matter 

We examine the ejection process of neutron-rich mat¬ 
ter in detail, which corresponds to the evolution of Ij,, 
focusing on representative ejecta (tracer particles). The 
time evolution of the Ij. for the selected tracer parti¬ 
cles is shown in Figure!^ based on B12/31.00 (prompt) 
and B11/30.25 (delayecfh respectively. All panels have 
the evolution history of Ij, for the two tracer particles 
labeled by T/,,nse and computed by TPM. Each particle 
has two different lines, where the thick line shows evolu¬ 
tion of the Yf. including all weak interactions described in 
Section [2.3.2[ and the thin line indicates the result of the 
case ignoring neutrino absorption. We plot the cases of 
prompt-jet explosions in (a) and (b) of Figure^ These 
tracer particles correspond to very neutron-ricn matter 
(a). Ye < 0.2, as well as mild neutron-rich (b), Fj. ~ 0.3, 


respectively. The evolution of Ye based on the delayed 
explosion model is also plotted in (c), where the final 
Fe > 0.3. 

We can qualitatively understand the role of neutrino 
absorption on the evolution of F),. For all the cases, the Ye 
of the ejecta rises via the neutrino absorption, where the 
reaction n -|- t'e —t p -I- e“ is predominant for neutron-rich 
(Fe < 0.5). The activity of this reaction is determined by 
the distance from the PNS (neutrino-sphere), and the lu¬ 
minosity and mean energy of the emitted neutrino, which 
appear to be independent of the geometry of explosion 
(the direction of ejection). Thus, during the expansion 
phase, the effect of neutrino absorption is well explained, 
assuming the spherical symmetry geometry, where reac¬ 
tion rates are just determined by the distance, even in 
cases of multi-dimensional hydrodynamics. Each tracer 
particle increases the Fe (thick line) due to neutrino ab¬ 
sorption up to three times compared with the case with¬ 
out neutrino absorption (thin line). The physical uncer¬ 
tainty of neutrino reactions due to the properties of PNS 
is discussed in the following section. 

The neutron-richness, which is in inverse proportion 
to the Ye, is determined by the strength of the electron 
captures on protons, which become stronger in higher 
densities. In order to examine more details, we now de¬ 
scribe the hydrodynamical evolution of the temperature, 
density, and radial distance from the center for the tracer 
particles shown in Figure Neutron-rich ejecta in the 
prompt-jet model, Fg^nse = 0.152 and 0.198 in Figure]^ 
and left panel of Figure ^ have higher maximum densi¬ 
ties (> 10^^ g cm^) at the bounce}^ where the electron 
capture is more active than other weak interactions. On 
the other hand, the particles, Fg^nse = 0.295 and 0.315 in 
Figure [^, have lower densities in the early stage of the 
expansion. These F/,s increase after the bounce, which 
implies that the effect of neutrino absorption is stronger 
than electron captures. 

The behavior of the Fj. evolution for the delayed model 
appears to be rather complicated. The ejecta initially 
stays in the PNS core, which has a higher density (~ 
10^^ g cm“^) and temperature (~ 10^^ K), whose Fg is 
determined by the /3-equilibrium condition. These parti¬ 
cles remain in the core until 0.05 s (= 50 ms) after the 
(failed) bounce, where this time duration corresponds to 
the timescale of the magnetic fields enhancement (tdei) 
in Table 15 The density and temperature immediately 
drop, resmting in the lower densities, during the later 
phase of the explosion. The Fg sharply increases at the 
beginning of the expansion, due to strong neutrino ab¬ 
sorption, beca use they are close to t he neutrino-sphere 
(also found by |Winteler et al.||2012|). Here, we should 
note that the effects of neutrino absorption (weak in¬ 
teraction) are sensitive to the treatment of neutrino re¬ 
actions and the transport, especially for the delayed-jet 
models. More sophisticated treatments are desirable for 
more precise predictions for the Fj. of ejecta. 

The maximum density and temperature of the tracer 
particles with the F(,_nse are plotted in Figure for 
Bll/31.00 (prompt) and Bll/30.25 (delayed), respec¬ 
tively. The time of the maximum density almost cor¬ 
responds the time of the maximum temperature, which 


In this context, we use the term “bounce” for each particle, 
which does not correspond to the time of bounce at the PNS core. 
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Fig. 10.— The histogram of ye,nse vs. the mass fraction of the total ejecta. Each plot corresponds to the explosion model whose name 
is labeled in the upper left, where there are two different histograms, that is, the case includes neutrino absorptions (red bars) and ignores 
these effects (blue bars). 
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Fig. 11.— The histogram of Ye for all ejecta changin g th e radius of the neutrino-sphere where the mass fraction is plotted in linear scale. 
Each panel has two bars in addition to those in Figure |l0| (green for the delayed model and red for prompt models), for which we assume 
a longer (blue) and shorter (magenta) time duration of neutrino absorptions, respectively. 
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is also close to the bounce time of each particle. In each 
panel of Figure particles in higher densities and tem¬ 
peratures (upper right) correspond to the ejecta that has 
escaped from deepest regions inside the PNS, of which 
the matter finally remains neutron-rich (T(j^nse < 0.3). A 
group of particles in the low densities and temperatures 
(lower left), on the other hand, is composed of ejecta in 
the jet-cocoon pushed by the jet-like outflow without col¬ 
lapse. The latter component of the ejecta remains initial 
Ye (> 0.4) at the pre-collapse stage, where the r-process 
hardly occurs. 

The values of Yg for neutron-rich matter are differ¬ 
ent between the prompt and delayed models, although 
they have a similar distribution of the maximum den¬ 
sity and temperature (Figure]^. As we have described, 
prompt-jets have very neutron-rich ejecta of Fe,nse < 0.2, 
whereas the ejecta of the delayed-jet hardly falls below 
k)i,nse 0.3. This is caused by the difference of the ejec¬ 
tion processes after the core-bounce, which is indepen¬ 
dent from the maximum values of densities and temper¬ 
atures. Here, we recall that the timescale for the amplifi¬ 
cation of magnetic fields (tdei in Table or the duration 
of magnetic fields wrapping) is a crucim factor to deter¬ 
mine the property of the explosion models around the 
PNS. Additionally, the delayed model has a lower explo¬ 
sion energy, in which the densities and temperatures are 
lower, leading to weaker electron captures on the pro¬ 
tons. The delayed model in the present study, which 
is based on a long-term SR-MHD simulation, has been 
overlooked in previous nucleosynthesis studies because 
they mostly focused on strong magnetic driven jets cor¬ 
responding to prompt-jets in the present classification. 
However, delayed-jet models, which are expected to have 
different elemental products, cannot be ignored for the 
nucleosynthesis study of MR-SNe. Finally, we should 
note that the property of the PNS is sensitive to relevant 
physical processes. We evaluate the uncertainty of Yg 
based on our current explosion models in the following 
section. 


3.2.3. The Yb of ejecta and uncertainty of PNSs 

We have calculated the evolution of Yg using the TPM 
for all our explosion models. Figure [l0| shows the his¬ 
tograms of the Yg (at the end of NSE) for the entire 
ejecta in a range of 0 < Fg, nse < 0.45, where the width 
of Yg (Fg-bin) corresponds to AFg = 0.02. Each model 
has two different Tg-histograms, which take into account 
all relevant weak reactions (green or red bars of the upper 
panel) and ignore neutrino absorption during expansion 
(blue bars of the lower panel), respectively. 

Prompt-jets (red histograms in Figure |10[ ) have lower 
Yg ejecta compared with the delayed mooel (green), of 
which peak at T(j,nse 0.3. The prompt-jet has a variety 
of Yg distribution among all explosion models. We can 
confirm that all prompt models have very neutron-rich 
ejecta of Fg^nse < 0.2, and expect successful r-processes 
producing the third peak. As we show in the following 
nucleosynthesis calculations, the difference in the Yg^nse 
among all prompt-jet models becomes small when we see 
the resulting nucleosynthesis abundance patterns. On 
the other hand, the delayed model (the green histogram 
in Figure 10) mostly has moderate neutron-rich ejecta 
Fe,nse ~ O.ST^of which the production of heavy r-process 
elements is restricted. This difference of the Fg-histogram 


between the prompt and delayed models remains signif¬ 
icant in the final abundance pattern. 

Here, we will emphasise the importance of neutrino ab¬ 
sorption in MR-SN models for determining of Fg. In Fig¬ 
ure the histogram in the lower panel of each model 
shows result that we ignored the effect of neutrino ab¬ 
sorption during expansion (the neutrino-optically thin 
region). The Fg generally increases due to neutrino ab¬ 
sorption in realistic explosion conditions, in which we 
have shown the time evolution of selected tracer particles 
in Figure By comparing these two Ij^-histograms, we 
find that the peak of extremely low 1), < 0.1 in the case 
of ignoring neutrino absorption is removed or shifted to 
higher Fg (the difference is AFg ^ 0.1) when we take 
into account neutrino absorptions. Although detailed 
treatments of weak reactions are essential to calculate Fg, 
the precise treatments of neutrino transport fully coupled 
with multi-dimension hydrodynamics are still computa¬ 
tionally expensive. 

Therefore, we take into account the physical uncer¬ 
tainties relevant to the weak interaction and the PNS 
by changing the radius of the neutrino-sphere when we 
calculate the evolution of Yg in the TPM. Based on the 
standard value, we assume a larger and smaller radius, 
which corresponds to the shorter and longer time dura¬ 
tion of neutrino absorption, respectively. The fiducial 
values of the radius are generally ^ 100, as shown 
in Table 1^ where we assume the variation of ±20% (i.e., 
the differences are ~ 10 km). 

The result of Fg,nse, based on the the uncertainty of the 
PNS, is shown in Figure [Tl] which plots the histogram 
of Fg for all explosion models. They generally exhibit a 
lower Yg distribution for a smaller radius case, whereas a 
larger radius of the neutrino-sphere results in higher val¬ 
ues for the Yg. This is because a larger neutrino-sphere 
causes ea^ ejection from the neutrino-sphere (Phase 1 
in Figure^, which means electron captures reducing the 
Yg are active. The variation of the neutrino-sphere ra¬ 
dius also affects the reaction rates of the neutrino inter¬ 
action outside the neutrino-sphere, in which the smaller 
radius results in weaker neutrino absorption. In contrast, 
the case of the larger neutrino-sphere radius has the 
longer duration staying the neutrino-sphere (i.e., tracer 
particles are influenced by weaker electron captures and 
stronger neutrino absorption). The delayed model is ba¬ 
sically affected by neutrino absorption, so that the effect 
of uncertainty o n t he Fg is larger than prompt-jets, as 
shown in Figure [TT| 

4. NUCLEOSYNTHESIS 

Nucleosynthesis calculations are performed based on 
the MR-SN models described in the previous section. 
The results are compared with several r-process abun¬ 
dance patterns of the sun and metal-poor stars. We dis¬ 
cuss the effects of physical uncertainties on nucleosynthe¬ 
sis, focusing on nuclear reactions, properties of PNS, and 
hydrodynamical instabilities of jets. Finally, the produc¬ 
tion of ®®Ni in jet-like explosions of the MR-SNe is briefly 
discussed with optical observation of SNe. 


4.1. The r-process Nucleosynthesis 

We performed nucleosynthesis calculations employing 
the nucl ear reaction network code described in Sec¬ 
tion 2.5 The network calculations are carried out by 
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Fig. 12.— The final abundances (Y) as a function of mass number (A) with the solar r-process pattern (open dots). Each explosion 
model has three different lines, which are the results of low (green), medium (blue), and high (red) resolution for adopted tracer particles. 


following the thermodynamical evolution of trajectories 
in the basis of the MR-SNe models. All of results for 
the explosion models are shown in Figure |12[ which 
are plotted integrated final abundances Y (the number 
abundance relative to the Si atoms, T(Si) = 10® ) as a 
function of mass number A of each isobar. Each panel 
has three lines for different resolution, which are labeled 
“low,” “medium,” and “high,” corresponding to the sets 
of tracer particles in Table Observational-based so¬ 
lar r-proces s abundances (the re sidual of the classical 
s-process by Arlandini et al.|[l9M I are depicted by open 
dots. 

We confirm that the nucleosynthesis calculations are 
numerically converged because the abundances of high 
and intermediate (medium) resolution are almost iden¬ 
tical. In the case of low resolution, we can see that 
some models are not converged. The models of B12/30.25 
and B12/I1.00, in particular, have a significant difference 
in abundance patterns between the low resolution and 
higher resolution calculations. In general, the number of 
tracer particles necessary to resolve the final r-process 
abundances depends on the hydrodynamical properties 
of explosion models and the treatment of tracer parti¬ 
cle motion. Our results point out that several hundred 
particles (the current medium resolution) in 2D hydrody¬ 
namics simulations are required to resolve the resulting 
r-process, which has not been satisfied in previous stud¬ 
ies. 

The mass of r-process elements in the ejecta and the 
total ejected masses of a jet are summarized in Table 
for all explosion models. The total mass of the ejecta is 
estimated only for the jet-like outflow, which is located 
in < 4000 km at the pre-collapse phase. The total mass 
of the r-process material is in the range of 0.957 x 10“^ 
- 2.45 X 1 O“^M0, where the less energetic delayed model 
of Bll/30.25 shows a relatively lower value and the most 


TABLE 5 

Masses of r-PROCESS Elements and the Total Ejecta 


Model 

r-process Mass 

(10-2m«) 

Ejected Mass 
(lO-^Mn) 

Bll/30.25 

0.957 

2.68 

Bll/31.00 

1.53 

2.15 

B12/30.25 

1.26 

3.55 

B12/31.00 

1.66 

4.37 

B12/34.00 

2.45 

8.57 


Note. — The amount of r-process elements in the ejecta and 
the total ejected mass in the solar mass for our MR-SN explosion 
models. 


energetic prompt model, B12/34.00, has a much larger 
value. These higher masses are necessary for ejecting 
very neutron-rich matter from the deeper region of a 
PNS, avoiding the Vg equilibrium, which is different from 
the moderate neutron-rich ejecta of PNS winds. Using 
an analytic formula of the minimu m mass loss rate for a 
neutron-rich ejecta (Equation 19 of [Metzger et al.|2008 ), 
we can derive ^ 1 x 10 “^Mq s“^ adopting typical values 
of Tablewhere ~ 4 x 10®^ erg s“^, ~ 80 km, 

and ~ 15 MeV with the 1.4Mq PNS mass. This mass 
loss rate results in ^ 1 x 10 “^Mq, with the duration of 
mass ejection of 200 ms, which is consistent with results 
from the hydrodynamics simulations. 

The higher mass ejection of the r-process elements sup¬ 
ports the importance of MR-SNe in the galactic chemical 
evolution. The typical ejected mass of r-process nuclei 
in MR-SNe, ~ 1O“^M0, is several thousand times larger 
than the typical ejected mass of PNS winds (the case of 
canonical CC-SNe). These higher values of r-process el¬ 
ements indicate that MR-SNe have a significant impact 
on the galactic chemical evolution if their event rate is 
very low. If we assume that only MR-SNe are the source 
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Fig. 13.— The final abundances of a ll explosion models ar e 
plotted with the solar r-process pattern l|Arlandini et al.||1999||. 
Red lines indicate the results for the prompt-jet and ttie green line 
corresponds to the delayed-jet. 

of r-process elements, we can estimate that ~ 0.1% of 
all CC-SNe explain the amount of the Galactic r-process 
material by simple multiplication. 

4.1.1. The Prompt-jet versus the Delayed-jet 

Comparing th e fin al abundance patterns of all models 
shown in Figure [T^ we find common global features in 
prompt-jets; these are all models with the exception of 
Bll/30.25, which is the only delayed model. Prompt- 
jet models, which are Bll/31.00, B12,50.25, B12/31.00, 
and B12/34.00, successfully produce heavy r-process nu¬ 
clei, including the second and third peaks and actinides, 
although the abundances of lighter isotopes have disper¬ 
sion. On the other hand, the final abundance curve of the 
delayed model (Bll/30.25) reaches up to the second r- 
process peak, where heavier nuclei {A > 100) are severely 
underproduced. These similarities and the difference in 
the final abundances among our M R-SN m odels are con¬ 
sistent with the discussion in Section 3.2.3 which is based 
on the Te-histograms in Figure [TO] 

The final abundances of Figure |13| are the results of the 
nucleosynthesis calculations base3~on the same nuclear 
physics input (i.e., experimental and theoretical masses, 
reaction rates, and decay properties). Therefore, each 
nucleosynthesis abundance pattern reflects the property 
of each hydrodynamical model, determined by stellar ro¬ 
tation and magnetic fields. Our results clearly indicate 
that the rotation and magnetic fields of stellar models 
change the r-process nucleosynthesis, that is, the nucle¬ 
osynthesis signature of the MR-SNe has a wide variety 
of observable features (chemical abundances), caused by 
the dynamical properties of stellar evolution models. In 
particular, the prompt and delayed models appear to 
be distinguished in the characteristics of the r-process 
yields. 

All of the prompt-jet models, which are energetic jet¬ 
like explosions, have the similar feature of the r-process 
yields, although the initial conditions and the details of 
the following explosion process are different. On the 
other hand, the weaker explosion of the delayed-jet pro¬ 
duces elements up to the second peak. The physical 
properties of MR-SN explosion models continuously vary, 
depending on the initial rotation and magnetic fields. 


Thus, the transition of nucleosynthesis yields from the 
prompt-jet to the delayed-jet appears to indicate the ex¬ 
istence of a threshold for the production of r-process 
nuclei. Although we found a qualitative difference be¬ 
tween the prompt and delayed-jets, the number of explo¬ 
sion models in the current study is insufficient to identify 
these thresholds. Further systematic parameter studies 
are necessary for the clarification. 

4.1.2. Physical uncertainties 

The calculated r-process yields have uncertainty due 
to the nuclear physics input for the reaction network. As 
described in the previous section, the prompt-jet repro¬ 
duces solar r-process abundances well, while there there 
exists deviations around the second and the third peak 
and an imderproduction of the rare-earth peak region 
(Figure [T^. It has been pointed out that the deficiency 
is mostly cause d by nuclear physics uncertainty (com¬ 
pare with, e.g., [Winteler et al.||20I2| utilizing the same 
FRDM (1995) mass model). Impro vement of this mass 


model, which is a new FRDM 
[2M4|, 

rROT 


(see, 


Kratz et al. 


Goriefy et a 


Moller et al. 


or di fferent mass models (e.g. 


2012 

HFE 


pDU9 |2010|| are expected to solve this prob- 


lem. Resolving these nuclear physics uncertainties is im¬ 
portant, however, as long as we focus on global features 
of abundance patterns (e.g., the relative amounts of the 
second and third peaks and the heaviest nuclei produced 
r-processes), the difference due to explosion dynamics 
is much larger. Therefore, the current discussion about 
the r-process in the MR-SN models (i.e., the impact of 
explosion dynamics on nucleosynthesis), is expected to 
be applicable to the results based on different reaction 
networks, as demonstrated in previous stud ies (see, e.g., 
Nishimura et al.||2006 Fujimoto et al.||2008). 

Aside from the nuclear physics inputs of the nuclear 
reaction network, the physics of PNSs have additional 
uncertainties. The evolution of neutron-richness (or Fj.) 
during the core-collapse, the core-bounce, and an early 
phase of ejection is i nfluenced by weak interactions. As 
shown in Figure 0 the Ye of ejecta has a significant 
change within a reasonable range of variation for the ra¬ 
dius of a PNS, adopting ±20% of the standard value. 
Following this evaluation, we recalculate r-process nucle- 


osynthesis, of which the results are shown in Figure 14 
Each explosion model has a final abundance of the “stan¬ 


dard Ye” (the case of high resolution in Figure 12), with 
different sets of Fe’s that corre spond to the smalFT?^ and 
large-(described in Section |3. 2. 3|), respectively. 

For prompt-jet models (Bll^l .00, B12/30.25, 

B12/31.00, and B12/34.00), the global feature of the 
abundance patterns remains, although the effects of 
the Ye variation on the final abundances appear in 
the amount of the rare-Earth peak, actinides, and 
digs around the second and third peaks. We should 
recall that these isotopes are also sensitive to nuclear 
physics inputs for the reaction network calculations. 
Therefore, we can argue that the prompt-jet of MR-SNe 
produces and ejects heavy r-process nuclei, reproducing 
the solar-abundance curve within the uncertainty of Yg. 
However, in order to determine the more precise value 
of Ye and the final abundances, we need further studies, 
based on more sophisticated treatment of the neutrino 
transport. 

On the other hand, the delayed-jet Bll/30.25 has larger 
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Fig. 14.— The final abundances (Y) as a function of mass number (A), based on the different set of Ye, wi th t he solar r-process pattern 
(open dots). Each explosion model has three different lines, following distribution of Ye shown in Figure [Tu which are the results of 
standard (green for the delayed model or red for the prompt model), small-(blue), and large-iljy (magentaj, respectively. 


effects on the uncertainty of the PNS on the final prod¬ 
ucts. This is because of a larger variation of com¬ 
pared with prompt models due to neutrino absorption. 
This is consistent with our discu ssion on the progress 
of the r-process in Section 3.2.3 only focusing on the 
neutron-richness of tracer particles. The abundance pat¬ 
tern of the delayed model reaches up to A ~ 130 for the 
standard case, so that the qualitative difference in the 
region of nuclei A > 130 appears to be significant. The 
case of small-of which the average value of Ye is low, 
slightly produces heavier r-process nuclei {A > 120), al¬ 
though this still underproduces compared with the solar 
abundances. This difference becomes important when we 
compare the final abundance with weak r-process pat¬ 
terns in the following section. 

4.1.3. The effect of Jet-deformation 

MR-SNe models employed in this study are based on 
2D axisymmetric simulations, which ignore several hy- 
drodynamical instabili ties that possibly o ccur in 3D. In 
particular, as shown by Mosta et al. (2014|), the launched 
jet-like can be deformed and destroyed via a hvdrodv- 
namical instability (e.g., the kink instability) The 
turbulent motion of deformed jets can induce reconnec¬ 
tion of the strong magnetic fields and part of the mag¬ 
netic energy can be converted to thermal energy. This 
heating raises the temperature during the expansion and 
can possibly destroy heavy elements that were made in 
the early jet ejection. An increase in the initial entropy 
for the same Yg is not an obstacle for r-process produc- 
tion, in fact, it provides m ore promising conditions (see, 
Freiburghaus et al.|1999a|. But, if such products experi¬ 


ence a strong heating after their initial production, these 
heavy elements can be destroyed by photo-disintegration. 

Here, we examine the effect of jet deformation on the 
nucleosynthesis, with a simplified evaluation of the re¬ 
heating due to reconnection. We select two different 
represented trajectories from the prompt-jet model of 
B11/31.00, whose temperature evolutions are shown in 
Figure [l^ The corresponding Ye at the end of the NSE 
phase are 0.180 and 0.326, leading to the production 
of heavy r-process nuclei, including the third peak plus 
lighter elements up to A ^ 130, respectively. For the lo¬ 
cation of the energy injection by magnetic reconnection, 
we chose three different cases at 300, 500, and 800 km 
from the c enter with typical mag netic field strengths (for 
details, see Takiwaki et al.|2009 ) of 1.0 x 10^"^, 5.0 x 10^^, 
and 2.0 x 10^'^ G, respectively. At the moment of the 
deformation, we assume that 25 or 50 % of the total 
magnetic energy Amag is converted instantaneously to 
thermal energy, as denoted by Em- Based on the in¬ 
creased total thermal energy, we estimated the increase 
in temperature using the Timmes EoSp^For the cases of 
larger energy deposition (e.g., Ain = O.bAmag at 800 km), 
the entropies are increased by 4.09 and 18.5 fee baryon” 
for F(j_nse = 0.180 and F(j_nse = 0.326 trajectories, respec¬ 
tively. 

The results of nucleosynthesis calculations using the 
modified trajectories are shown in Figure 16 For the 


This, of course, also depends on the exact rotation rate and 
magnetic fie ld strength, see also the 3D simulations of |Winteler| 


|et ai.lpM^. 


cases of Ye — 0.180 with larger energy deposition Ain = 
0.5Aniag (a), the production of heavier nuclei is partially 
suppressed. In particular, the model with energy injec¬ 
tion at 800 km shows the largest difference, changing 
entire r-process nuclei by a factor, although the main fea¬ 
tures of the pattern still remain. For the lower injection 

http://cococubed.asu.edu 
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Fig. 15.— The evolution of temperature in GK for selected tracer 
particles, labeled by the value of Ye at the end of NSE. Heating 
due to the magnetic reconnection is taken into account at 300, 500, 
and 800 km from the center, respectively. 

energy Ei^ = 0.25£'niag (b), all models show a smaller 
impact by the heating, except for the isotopes A > 200. 
For the Ye = 0.326 trajectory with large Ein = 0.5£'niag, 
which produces mostly lighter r-process nuclei A < 130, 
most of the isotopes with A < 100 are not changed. Sim¬ 
ilar to case (a) above, the heavier nuclei with A > 100 are 
somewhat suppressed by the late temperature increase. 

Based on the estimation of the simplified models, we 
can conclude that in the case of very strong energy de¬ 
position due to magnetic reconnection (e.g., if kink in¬ 
stabilities occur), the late temperature increase can lead 
to the destruction of heavy elements after their initial 
production. Therefore, the suppression of the heaviest 
elements can take place (reducing them by a factor), but 
these effects do not change the main feature of our cur¬ 
rent results. 

4.2. r-process Enhanced Metal-poor stars 

In addition to the solar abundances, we compare our 
results with the observed abundances of metal-poor stars. 
In Figure results of the calculations are shown with 
the abundance patterns of two r-process-rich metal-poor 
stars. We choose Bll/31.00 as a typical prompt-jet 
model, while we adopt the standard set of Ye and the case 
of small-for the delayed-jet model Bll/30.25. These 


final abu nda nce patterns are based on the same result 
in Figure [131 although we replotted as a function of the 
atomic number Z. In this figure, these are compared 
with the abundance patterns of two different r-process- 
enha nced metal-poor star s, which we adopt CS22892-052 


(e.g., Sneden et al. 19961 and HD122563 (Honda et al. 
2004bF I. They represent solar-like r-process abundances 


and the “weak r-process” pattern, respectively^^ The 
solar-like distribution has been observed in severaTmetal- 
poor stars, which have a remarkably similar feature, or 
“university,” for heavier r-process elements {A > 130). 
So far there have not been as many examples repo rted 
as weak r-process patterns (see, Sneden et al.||200^ . 

The final abundances of pronipt-jet models are basi¬ 
cally agree with the observational pattern of CS22892- 
052, which is a solar-like r-process abundance distribu¬ 
tion, although the underproduction is around Z = 60. 
This is of course same as the case of the comparison with 
the solar-abundance pattern, as discussed previously. In 
the case of the delayed-jet, on the other hand, the result 
of standard Ye produces nuclei up to Z ~ 55 and has 
no production for further heavy elements. As we con¬ 
sider the uncertainty of 1^, this severe underproduction 
appears to be dissolved, because the delayed model with 
smaller R^, (smaller Ye) produces heavier r-process nu¬ 
clei. The r-process in the low entropy environment is 
much depend on the degree of neutron-richness, so that 
further lower Yg cases, which are between the prompt 
and delayed-jet models, are expected to be reproduce 
the entire pattern of HD122563. 

When we consider MR-SNe as the origin of nucleosyn- 
thetic signatures in r-process-rich metal-poor stars, we 
can repeat the discussion on the solar abundances. The 
ejecta of prompt-jet models can be a source for CS22892- 
052, which is a solar-like r-process-rich metal-poor star. 
Whereas, the weak r-process pattern seems to indicate 
an intermediate feature between our prompt and delayed 
models. The question then arises whether its features re¬ 
late to a single/pure process from one event or whether 
it is a wei ghted superposition of different events (as dis¬ 
cussed by Hansen et al.|2014 ). In future, more extended 
parameter studies, we will explore whether the full weak 
r-process pattern can be reproduced in total. This might 
offer the chance to explain the scatter and different r- 
process abundance patterns in the early galaxies. 


4.3. synthesis oj Ni isotopes in Jets 

In addition to r-process nuclei, CC-SNe produce a large 
amount of iron-group isotopes by explosive nucleosynthe¬ 
sis. Radioactive nickel isotopes (e.g., ®®Ni) are of partic¬ 
ular importance for connecting nucleos ynthesis yields of 
SNe t o several optical observations (see,|Magkotsios et ah’ 
2010 and references therein). For instance, emission of 
7 -rays from ®®Ni via the chain ^®Ni — >• ®®Co ®®Fe is 
the dominant energy source of the early phase of SN light 
curves. The jet-like explosion of MR-SNe ejects not only 
very neutron-rich material from the inner region, but also 
outer layers of the iron core Yg > 0.4 as shown in Fig¬ 
ure 1^ Here, we calculated the ejected masses of ®®Ni, 


There are still several arguments of whether or not a single 
event originates this pattern. In this paper, we mention it as an 
alternative observational r-process signature. 
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Fig. 16.— Final abundances of different energy deposition location of the magnetic reconnection at 300, 500, and 800 km from the 
center: (a) the trajectory of ye,nse = 0.180 with the energy injection E[n = O.S^mag- (b) Same as (a), but the lower injection energy 
E[n = O.^S-Emag- (c) The trajectory of y’e,nse = 0.326 with the injection energy of Ein = O.SEmag- 



Fig. 17.— Final abundances as a function of atomic number Z 
with the abundances of r-process-rich metal-poor stars of CS22892- 
052 and HD122563. The solid line (red), and the thick dashed 
line (green), and the thin dashed line (green) are a prompt-jet 
model (B11/91.00) and a delayed-jet model (Bll^O.25) and a de¬ 
layed model with small-respectively. 


®^Ni, and in the primary MHD-driven jet which 
is summarized in Table Synthesis of these isotopes 
takes place in a distant region from the PNS, where the 
effect of neutrino absorption is negligible. In fact, the 
final abundances of iron-group elements change only a 
few per cent, b ased on the uncertainty of Ij, discussed in 
Section 14.1.21 

The mass of ®®Ni, M(®®Ni), is in the range of 0.714 x 
1 O“^M0 and 4.27 x Mq, which appears to be cor¬ 
related with the total ejected mass in Table Our 
models have lower values compared with a typical type 
SN Up, SN 1987A, = 7.1 x and 


M(°' Ni) = 4.1 X 10 ^Mq (see, e.g., Seitenzahl et al. 


20141. The jet-like explosion of MR-SNe sometimes has 


been considered to be the central engine of luminous SNe 
Ib/c or hypernovae associated with long GRBs. How¬ 
ever, this amount of ®®Ni can hardly explain these lumi¬ 
nous events, which require M(®®Ni) > O.IMq. While, 
M(®®Ni) ^ O.OIMq corresponds to the lower bound of 
ordinary SNe, which is in the faint SN branch (see, Fig- 


We excluded B12/90.25, which does not yet feature the full 
explosive nucleosynthesis. The current explosion models have been 
geared mainly to investigate the r-process, as determined by the 
Ye at the early phase of explosions. 


TABLE 6 

Masses of nickel isotopes 


Model 

M(5«Ni) 

(lO-^Me) 

M(5^Ni) 

(lO-^Me) 

M(5*Ni) 

(10-3) 

Bll/30.25 

1.08 

2.51 

0.774 

B12/30.25 

0.714 

2.48 

3.31 

B12/31.00 

1.06 

2.91 

2.57 

B12/34.00 

4.27 

7.81 

6.02 


Note. — The amounts of ^®Ni, ^^Ni, and ^®Ni in the jet¬ 
like ejecta, as measured in the solar mass (Mq). The model of 
B12/30.25 is excluded because the duration of the simulation is in¬ 
sufficient for determining these isotopes. 


ure 1 of Nomoto et al.|[2006 Smart! 
' 2009[ ). This indicates that the de' 


2009 Tanaka et al.l 


ayed-jet model, of 
which ^®Ni mass is ~ O.OIMq with a less energetic ex¬ 
plosion, may be classified in the faint SN. 

Based on the amounts of Ni isotopes, we can quan¬ 
titatively discuss the connection between MR-SNe and 
several high-energy astronomical phenomena. An X- 
ray flash event, XRF 060218, that is associated with 
SN 2006aj (Type Ic) requires ®®Ni ~ 0.2Mp) and ^ ^Ni 


0.05A/f7^ by optical observations (Mazz ali et al. 2006 
Maeda et al.||2007a| Mazzali et al. 20071. AlthouA oui' 
explosion models show lower INi-isotope ejecta, they are 
still suitable for another peculiar type Ib SN, SN 2005bf, 
which has intermediate explosion energy between normal 
type Ib/c SNe and X-ray flashes. Our results have the 
amount of Ni isotopes in the range of M (®®Ni) < O.OSMq 
with ~ 10^^ - 10^® G, which are values for a magnetar 
formation scenario (Maeda et al. |2007b| . Additionally, 
the lower ejected mass of M(^'^Ni) = 0.00 3Mf7^ is sug¬ 


geste d for jet-like explosion of SN 2010jp (Smith et al. 
|2012| ), which may correspond to a less energetic MR-SiN 
mo dels 

However, luminous SNe Ib/c and relevant explosive 
phenomena including GRBs and XRFs gen erally requir e 
larger amounts of ®®Ni > 0.2 Mq (see, e.g., Gano|[20T3 ), 
which the calculated value of ®®Ni for our explosion mod¬ 
els hardly explain. Here, it should be emphasized that 
we have only focused on the production of Ni isotopes in 


The collapse model also may be suitable for such low 
ejecta (the lower bound of ®®Ni is 3.7 X by|Fujimoto et al.| 

|2008||. However, it is not clear if BH-driven jets are possible with¬ 
out polluti on of yields by magn etar-driven (MR-SN) jets before BH 
formation ([Harikae et al.|2009|. 
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the primary jets. Thus, we have ignored additional ejecta 
due to neutr ino-heating after jet prop agation, which is 
discussed by |Sawai & Yamada| (|2014 ). We can expect 
significant arnounts ot JNi and other iron-group isotopes 
in this neutrino-driven ejecta, unless the central PNS im¬ 
mediately collapses to BH. The current study might give 
a clue for further investigation using more sophisticated 
models. 


5. SUMMARY AND CONCLUSIONS 

We investigated the r-process nucleosynthesis, focusing 
on the explosion scenario of MR-SNe. We found char¬ 
acteristic nucleosynthesis signatures, which are different 
from isotopic abundances in canonical neutrino-heating 
driven CC-SNe. The results are summarized as follows: 

1. A successful r-process, producing the solar r- 
process pattern, takes place in the jet-like ejecta 
of all prompt-magnetic-jet models, which have 
stronger initial magnetic fields. The final abun¬ 
dances are in good agr eement with solar r-process 
abundances (Figure fl^ within the remaining phys¬ 
ical uncertainties. 'The weaker magnetized explo¬ 
sion, categorized as the delayed-magnetic-jet, pro¬ 
duces matter up to the second peak nuclei {A ^ 
120 ). 

2. The prompt-jet can explain the solar r-process 
abundances and solar-like patterns of metal-poor 
stars, while the delayed-jets possibly contribute to 
“weak r-process” patterns of r-process-rich metal- 
poor stars. Additionally, our results predict that 
intermediate abundance patterns between a main 
(solar-like) r-process and a weak r-process are pro¬ 
duced by certain conditions of stellar rotation and 
magnetic fields. 

3. The total amount of ejecta enriched by r-process 
elements is in the range of ^ 10~^Mq per event. 
This value is about a thousand times larger than 
the typical value of PNS wind ejecta. Thus, MR- 
SNe have a significant contribution to the galactic 
chemical evolution of r-process elements, also re¬ 
quiring that the frequency of events is relatively 
low. This is consistent with the observational indi¬ 
cation of MR-SN events, of which frequency is low 
compared to canonical neutrino-driven SNe. 

4. The uncertainty of the present r-process calcula¬ 
tions due to the treatment of weak interactions is 
evaluated, focusing on the physics of the PNS. Rea¬ 
sonable uncertainties, which allow a ±20% varia¬ 
tion of the neutrino-sphere radius, modify nucle¬ 
osynthesis results, although the main properties 
found in this work remain. Delayed-jets, which are 
strongly influenced by neutrino absorption, experi¬ 
ence larger changes than prompt-jet models. 

5. We also estimated lower bounds for the amounts of 
Ni isotopes in jets. The amount of ^®Ni (~ O.OIMq 
for our explosion models) is significantly lower than 
the value in the ejecta of canonical CC-SNe (®®Ni ~ 
O.IM 0 ). Thus, primary jets of the MR-SN are ex¬ 
pected to be found as faint SNe by optical obser¬ 
vation. However, additional matter can be ejected 


in a later phase of the explosion, which can eject 
larger amounts of ®®Ni corresponding to luminous 
hypernovae and/or GRBs. 

The prompt-jet models exhibit similar nucleosynthesis 
features, in spite of different initial rotation and mag¬ 
netic fields. In other words, it is clear that strong polar 
jet-like explosions, in general, produce and eject signif¬ 
icant amounts of heavy r-process material. This indi¬ 
cates that the production of heavy elements is “satu¬ 
rated,” when initial rotation and magnetic fields exceed 
thresholds. Although our explosion models are based on 
the axisymmetric hydrodynamics, nucleosynthesis prop¬ 
erties of the prompt models are expected to be simi lar to 
the result of pola r-like jets in full 3D simulations (Win- 
teler et al. 2012). The impact of the remaining physi¬ 
cal uncertainties of the current explosion models on Yg 
(and the resulting r-process nucleosynthesis) are not neg¬ 
ligible, especially for the delayed-jet explosion, which is 
more sensitive to neutrino absorption. 

The amount of ejected r-process-rich matter exceeds 
^ 1 X 1 O“^M0, which is about a thousand times larger 
than masses of PNS winds. Therefore, the upper limit of 
the frequency of MR-SNe is estimated to be ~ 0.1% of 
all CC-SNe by simple multiplication. However, we need 
further observational restrictions to clarify the role of 
MR-SN events in the entire galactic chemical evolution. 
On the other hand, MR-SNe are more likely to be dom¬ 
inant sites of r -process nucleosynthesis for metal-poor 
stars (see, e.g., |Cescutti & Chiappini||2014t |Wehmeyeii 
et al.|[Ml5 ). In the early Calaxy, we expect MR-SMe to 
play the rriore significant role in comparison to other r- 
process events, such as compact binary mergers (before 
they set in), and it needs to be investigated whether their 
frequency is a function of metallicity. 

Additionally, MR-SNe also have an advantage to 
explain individual observed abundances of metal-poor 
stars , lik e CS22892-052 and HD122563 (discussed in Sec¬ 
tion 14/2). Our prompt-jet models reproduce solar-like 


patterns of metal poor-stars, while the delayed-jet model 
contributes a weak r-process pattern. The agreement of 
the heavy element production for prompt-jets of MR-SNe 
corresponds to the observational “universality” of heavy 
r-process elements in metal-poor stars. On the other 
hand, the “diversity” of r-process abundances in metal- 
poor stars for lighter r-process elements with A < 120 
appears to have different origins, and can be interpreted 
as the outcome of a variety of nucleosynthesis yields in 
MR-SNe due to different initial rotation and magnetic 
fields. Thus, we expect that MR-SN events contribute 
to the production of r-process nuclei more actively than 
other r-process sources in the early Galaxy. Additionally, 
the result that r-process yields of MR-SNe have signifi¬ 
cant variation due to the properties of progenitors (i.e., 
initial rotation and magnetic fields), is consistent with 
the need for multiple r-process astrophysical sites (dif¬ 
ferent abundance patterns), which have been suggested 
by at least two diffe rent process in resent investigations 
(Hansen et al.||2014). 

In addition to r-process elements, MR-SNe also pro¬ 
duce other radioactive nuclei (e.g., ^®Ni), which are 
sources of SN light curves. The amount of M(®®Ni) ~ 
O.OIM 0 in primary magnetic driven jets for our explo¬ 
sion models hardly explain the luminosity of SNe Ib/c 
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or hypernovae associated with long GRBs. This rela¬ 
tively small amount of ^®Ni corresponds to “faint SNe” or 
could explain magnetar birth events like SN 2005bf. Our 
explosion models eject insufficient amounts of ^®Ni and 
®®Ni to explain observations of SN 2006aj associated with 
XRF 060218. However, the current explosion models 
possibly underestimate these amounts, because we only 
take into account ejecta in jet-like explosions. Therefore, 
further comprehensive studies based on a wider range of 
parameters may show MR-SN models that satisfy these 
values. 

The explosion mechanism of MR-SNe has a remain¬ 
ing fundamental problem, which is the origin of such 
strong magnetic fields during the core-collaps e. One op¬ 
tion is the inheritan ce from stellar evolution (Maeder & 
Meynet 2004| 2005), which is a long-standing unsolved 
problernl The other option is the magneto-rotational 
instability (MRI), which is an enhancement process of 
magnetic fields in the central core during the collapse 
and explosion. The role of MRI during core-co llapse has 
been studied with a focus on several aspects ( Ahi^^anm 
et al. I [20031 lObergaulinger et al.||2009| [Masada et al. 


2012 lObergaulinger et al.|2014||Masada et al.|2015|). Re¬ 

cently, global MHU simulations with MKl reported suc¬ 
cessful jet-like explosions, starting fr om moderate mag¬ 
netic fields of pre-coll apse models (Sawai et al. 2013| 
Sawai fc Yamad^|2014 ). Although there are still several 
problems that should be solved (i.e., numerical conver¬ 
gence, 3D effects, and initial configurations of magnetic 
fields), the effects of the MRI are definitely important for 
further studies. 

A major aspect of the present investigations is the as¬ 
sumption of initial conditions with fast rotation rates and 
strong magnetic fields, which lead to an early and rapid 
jet formation and explosion. In reality it might be pos¬ 
sible that such strong magnetic fields do not exist at the 
onset, and they will only be generated in the PNS as a 
result of the MRI or a convective alpha-omega dynamo. 
This will happen on a longer timescale, leading initially 
to conditions like in a neutrino-driven explosion, which 
we then expect to be followed by a jet formation due to 
the amplified magnetic pressure. In such cases, matter is 
affected longer by neutrino radiation, which leads to an 
increase of Te, and the yields are expected to be similar to 
our “delayed-magnetic-jets,” unless the launched shock 


is violently strong enough to push out very neutron-rich 
matter from the deeper region of PNS. The answers to 
these open questions rely on improved stellar evolution 
calculations with clear predictions of rotation and mag¬ 
netic fields, as well as long-term MHD simulations in 3D, 
following the formation and evolution of the PNS. Such 
simulations are also expected to provide us with clues to 
connect the mechanism of MR-SNe with the shape and 
amplitude of the magnetic field in NSs. 

The nucleosynthesis properties have been presented in 
this paper show that MR-SNe play an important role in 
galactic chemical evolution and high-energy astronomical 
events. In particular, they may be the dominant source 
of r-process elements in early galaxies. The current hy- 
drodynamical simulations, however, are based on a sim¬ 
plified treatment of neutrino transport. Further studies 
with more sophisticated MR-SNe explosion models will 
hopefully provide more precise theoretical constraints 
and and predictions. Numerical data of trajectories and 
nucleosynthesis yields used in this study are available at 
http://www.astro.keele.ac.uk/~nobuya/mrsn. 
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